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EXECUTIVE  SUMMARY 

This  report  covers  work  accomplishted  at  the  University  of  Arizona,  Optical  Sciences 
Center,  during  the  period  from  December  15,  1993  through  May  14,  1997,  supported  by 
the  Joint  Services  Optics  Program.  This  final  report  describes  results  obtained  from 
eleven  research  projects,  spanning  the  field  of  optics  from  fimdamental  quantum  optics 
theory  to  the  engineering  design  of  optical  siystems.  These  projects  supported  more  than 
20  graduate  students  and  several  postdoctoral  associates,  and  the  research  effort 
culminated  in  the  publication  of  five  dissertations,  the  most  important  aspects  of  which 
are  listed  below.  They  have  been  published  by  University  Microfilms,  Inc.  (UMI),  and 
are  readily  available  from  UMI,  The  University  of  Arizona  Main  Library,  and  the  Optical 
Sciences  Center’s  Fred  A.  Hopf  Reading  Room. 

Thomas  C.  Zaugg,  “Cavity  QED:  Adiabatic  Atomic  Cooling  in  Cavities,  and 
Evaluation  of  a  Technique  for  Atomic  Homodyning  Detection  of  Cotangent  States,” 
University  of  Arizona  (1994).  This  thesis  analyzed  how  a  decaying  cavity  field  can  lead 
to  signficant  atomic  cooling.  It  also  includes  an  evaluation  of  a  nonlinear  atomic 
homodyne  detection  scheme  for  measuring  the  Wigner  characteistic  function  of  a 
microwave  cavity  field. 

Elena  V.  Goldstein,  “Nonlinear  Atom  Optics”,  University  of  Arizona  (1996).  The  goal 
of  this  resarch  was  to  study  selected  aspects  of  nonlinear  atom  optics.  Specifically,  the 
near-resonant  dipole-dipole  interaction  between  two  atoms  in  tailored  vacua;  the  study  of 
several  aspects  of  the  many-body  theory  of  atomic  ultracold  systems  in  situations  where 
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the  nonlinearity  arises  due  to  the  two-body  dipole-dipole  interaction;  and  the  elementary 
excitations  and  the  quasi-particle  frequency  spectrum  in  a  multicomponent  Bose 
condensate. 


J.  L.  Kann,  “Numerical  Modeling  of  a  Near-Field  Optical  System”,  University  of 
Arizona  (1995).  A  hybrid  finite-difference-time-domain  (FDTD)/angular  spectrum  code 
was  used  to  study  the  electromagnetic  and  imaging  properties  of  a  near-field  scanning 


optical  system. 

F.  F.  Froehlich,  “Resolution  and  Contrast  of  Optical  and  Shear  Force  Imagery  in  Near- 
Field  Scanning  Optical  Microscopy”,  University  of  Arizona  (1995).  Mechanisms  that 
influence  image  formation  in  near-field  scanning  optical  microscopy  performed  with 
tapered  fiber  aperture  probes  were  investigated. 


Byron  B.  Taylor,  “Topics  in  Atom  Optics”,  University  of  Arizona  (1997).  Investigations 
were  performed  on  the  effects  of  light  forces  on  the  center-of-mass  motion  of  two-level 
atoms,  Doppler  cooling  for  ry  optics,  and  an  atomic  Fabry-P6rot  for  wave  optics. 

The  projects  involved  contacts  and  collaborations  with  several  Department  of  Defense 


laboratories  including: 

U.S.  ARMY 

Natick  Resarch,  Development  and 
Enginering  Center 
Materials  Command 
Harry  Diamond  Laboratory 
Edgewood  Research,  Development 
and  Engineering  Center 


U.S.  AIRFORCE 

Phillips  Laboratory 
Rome  Air  Development  Center 
Brooks  Air  Force  Base 
Flight  Dynamics  Laboratory 
Wright  Patterson  Air  Force  Base 
Wright  Lab 


5 


A  complete  list  of  all  personnel,  publications,  and  inventions  are  listed  below; 

Personnel  Supported  by  the  Joint  Services  Optics  Program 

Galina  Khitrova,  Principal  Investigator,  “Nonlinear  Optics  in  Vapors  and  Electron 
Optics” 

Tom  Nelson,  Jr. 

David  Wick 
Eric  Lindmark 
Jill  Berger 
Ove  Lyngnes 

Hyatt  Gibbs/Michael  A.  Cusanovich,  Principal  Investigators,  “Photoactive  Yellow 
Protein” 

Gordon  Tollin 
Terry  Meyer 
Ove  Lyngnes 
Savitha  Devanathan 

C. F.  Li 

Pierre  Meystre,  Principal  Investigator,  “Atom  Optics” 

G.  Lenz 
T.  C.  Zaugg 
E.  V.  Goldstein 
M.G.  Moore 
P.  Pax 
B.  B.  Taylor 

Rudolph  H.  Binder,  Principal  Investigator,  “Theory  of  Ultrafast  Nonlinear  Optical 
Response  of  Uniaxially  Strained  Semiconductor  Quantum  Wells” 

Dariusz  Burak 
Nai  Kwong 
Gerald  Lapeyre 

Mahmoud  Fallahi,  Principal  Investigator,  “High  Power  Circular-Grating  Surface- 
Emitting  Lasers” 

Scott  Penner 

Poul  lessen.  Principal  Investigator,  “Atom  Trapping  in  Far-Off-Resonant  Optical 
Lattices” 

S.  E.  Hamann 

D.  L.  Haycock 
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G.  Klose 

E.  M.  Wright 
G.  J.  Lapeyere 
P.  H.  Pax 

Richard  Powell,  Principal  Investigator,  “Sodium  Guide  Star  Laser” 

James  Murray 

Ewan  Wright,  Principal  Investigator,  “Ultra-Short  Pulse  Propagation  and  the  Quantum 
Statistics  of  Optical  Solitons” 

Augusto  Rodrigues 
Keith  Kusanic 
Kirk  Cook 

Katherine  Creath,  Principal  Investigator,  “Real-Time  Electronic  Holographic  Moire” 
Norbert  Sigrist 
Joanna  Schmit 
Peter  Ettl 

John  Greivenkamp,  Principal  Investigator,  “Development  of  a  Generalized  Aspheric 
Interferometer” 

A.  E.  Lowman 

Tom  Milster,  Principal  Investigator,  “Micro-Optics  and  Near-Field  Scanning  Optical 
Microscopy  (NSOM)” 

J.  Kann 

F.  Froehlich 


GALINA  KHITROVA’S  PUBLICATIONS 

1.  H.  M.  Gibbs,  G.  Khitrova,  C.W.  Lowry,  and  S.  W.  Koch,  “Acceleration  of 
Coherent  Transfer  of  Energy  by  Stimulated  Emission  and  Absorption,”  Opt.  and 
Phot.  News  14  Dec.  1994). 

2.  R.  Jin,  G.  Khitrova,  D.  Boggavarapu,  H.  M.  Gibbs,  S.  W.  Koch,  M.  S.  Tobin,  and 
R.  P.  Leavitt,  “Physics  of  Semiconductor  Vertical-Cavity  Surface-Emitting 
Lasers,”  J.  Nonlin.  Opt.  Phys.  and  Matls.  4, 141-161  (1995).  Invited  paper. 


7 


3.  M.  F.  Krol,  S.  Ten,  B,  P.  McGinnis,  M.  J.  Hayduk,  G.  Khitrova,  and  N. 
Peyghambarian  “Subpicosecond  Hole  Tunneling  by  Nonresonant  Delocalization 
in  Asymmetric  Double  Quantum  Wells,”  Phys.  Rev.  B  Rapid  Communication  52 
(20),  14344  (1995). 

4.  P.  G.  Baranov,  N.  G.  Romanov,  I.  V.  Mashkov,  G.  Khitrova,  H.  M.  Gibbs,  and  0. 
Lyngnes,  “Local  Diagnostics  of  GaAs/AlAs  Superlattices  Using  Optically 
Detected  Magnetic  Resonance  and  level  Anticrossing  Effect,”  Solid  State  Physics 
37,1648(1995). 

5.  T.  R.  Nelson  Jr.,  E.  K.  Lindmark,  D.  V.  Wick,  K.  Tai,  G.  Khitrova,  and  H.  M. 
Gibbs  “Normal-Mode  Coupling  in  Planar  Semiconductor  Microcavities,”  in  J. 
Rarity  and  C.  Weisbuch,  eds.,  Quantum  Optics  in  Wavelength  Scale  Structures 
(NATO,  1995). 

6.  R.  Jin,  M.  S.  Tobin,  R.  P.  Leavitt,  H.  M.  Gibbs,  G.  Khitrova,  D.  Boggavarapu,  O. 
Lyngnes,  E.  Lindmark,  F.  Jahnke,  and  S.  W.  Koch,  “Order  of  Magnitude 
Enhanced  Spontaneous  Emission  from  Room-Temperature  Bulk  GaAs,”  in  J. 
Rarity  and  C.  Weisbuch,  eds..  Quantum  Optics  in  Wavelength  Scale  Structures 
(NATO,  1995). 

7.  G.  Mohs,  B.  Fluegel,  H.  Giessen,  K.  Meissner,  G.  Khitrova,  H.  Gibbs,  and  N. 
Peyghambarian,  “Ultrafast  Transient  Gain  in  Type  II  Multiple  Quantum  Wells,” 


Appl.  Phys.  Lett.  68  (4),  511  (1996). 
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8.  J.  Berger,  O.  Lyngnes,  H.  M.  Gibbs,  G.  Khitrova,  T.  R.  Nelson,  E.  K.  Lindmark, 
A.  V.  Kavokin,  M.  A.  Kaliteevski,  and  V.  V.  Zapasskii,  “Magnetic  Field 
Enhancement  of  the  Exciton-Polariton  Splitting  in  a  Semiconductor  Quantum- 
Well  Microcavity:  Beyond  the  Strong  Coupling  Threshold,”  Phys.  Rev.  B  54, 
1975  (1996) . 

9.  E.  K.  Lindmark,  T.  R.  Nelson,  Jr.,  G.  Khitrova,  H.  M.  Gibbs,  A.  V.  Kavokin,  and 
M.  A.  Kaliteevski,  “Three  Coupled  Oscillators:  Normal  Mode  Coupling  in  a 
Microcavity  with  Two  Different  Quantum  Wells,”  Opt.  Lett.  21, 994  (1996). 

10.  A.  V.  Kavokin,  M.  A.  Kaliteevski,  S.  V.  Goupalov,  J.  D.  Berger,  O.  Lyngnes,  H. 
M.  Gibbs,  G.  Khitrova,  A.  Ribayrol,  A.  Bellabchara,  P.  Lefebvre,  D.  Coquillat, 
and  J.  P.  Lascaray,  “Quantum  Wells  with  Zero  Valence  Band  Offset:  Drastic 
Enhancement  of  Forbidden  Excitonic  Transitions,”  Phys.  Rev.  B  54,  R1 1078 
(1996). 

11.  T.  R.  Nelson,  Jr.,  J.  P.  Prineas,  G.  Khitrova,  H.  M.  Gibbs,  J.  D.  Berger,  E.  K. 
Lindmark,  J.-H.  Shin,  H.-E.  Shin,  and  Y.-H.  Lee,  “Room-Temperature  Normal- 
Mode-Coupling  in  a  Semiconductor  Microcavity  Utilizing  Native-Oxide 
ALAs/GaAs  Mirrors,”  Appl.  Phys.  Lett.  69,  3031  (1996). 

12.  F.  Jahnke,  M.  Kira,  S.  W.  Koch,  G.  Khitrova,  E.  K.  Lindmark,  T.  R.  Nelson,  Jr., 
D.  V.  Wick,  J.  D.  Berger,  O.  Lyngnes,  H.  M.  Gibbs,  and  K.  Tai,  “Excitonic 
Nonlinearities  of  Semiconductor  Microcavities  in  the  Nonperturbative  Regime,” 


Phys.  Rev.  Lett.  77,  5257  (1996). 
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13.  S.  Ten,  J.  G  Williams,  P.  T.  Guerreiro,  G.  Khitrova,  and  N.  Peyghambarian, 
“Ultrafast  Excitonic  Room  Temperature  Nonlinearity  in  Neutron  Irradiated 
Quantum  Wells,”  AppLPhys.  Lett.  70,  158  (1997). 

14.  J.  D.  Berger,  S.  Hallstein,  0.  Lyngnes,  W.  W.  Ruble,  G.  Khitrova,  and  H.  M. 
Gibbs,  “Emission  Dynamics  of  a  Magnetoexciton  Quantum  Dot  Microcavity 
Laser,”  Phys.  Rev.  B  55,  R4910  (1997). 

15.  R.  I.  Dzhioev,  H.  M.  Gibbs,  E.  L.  Ivchenko,  G.  BChitrova,  V.  L.  Korenev,  M.  N. 
Tkachuk,  and  B.  P.  Zakharchenya,  “Linear-Circular  Polarization  Conversion 
Under  Optical  Orientation  of  Excitons  in  Type-II  GaAs/AlAs  Superlattices,” 
submitted  to  Phys.  Rev.  B. 

16.  M.  V.  Belousov,  A.  Yu.  Chernyshov,  I.  V.  Ignatev,  I.  E.  Kozin,  A.  V.  Kavokin,  H. 
M.  Gibbs,  and  G.  Khitrova,  “Statistical  Model  Explaining  the  Fine  Structure  and 
Interface  Preference  of  Localized  Excitons  in  Type-II  GaAs/AlAs  Superlattices,” 
submitted  to  Phys.  Rev.  B. 

17.  S.  Hallstein,  J.  D.  Berger,  M.  Hilpert,  H.  C.  Schneider,  W.  W.  Ruble,  F.  Jahnke,  S. 
W.  Koch,  H.  M.  Gibbs,  G.  Khitrova,  and  M.  Oestreich,  “Circularly  Polarized 
Laser  Emission  Locked  to  the  Electron  Spin  Precession  Clock,”  submitted  to 
Phys.  Rev.  Lett. 

18.  O.  Lyngnes,  J.  D.  Berger,  J.  P.  Prineas,  S.  Park,  G.  Khitrova,  H.  M.  Gibbs,  F. 
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HYATT  GIBBS’  PUBLICATIONS 

1 .  Lyngnes,  E.  K.  Lindmark,  M.  Switala,  H.  M.  Gibbs,  T.  E.  Meyer,  G.  Tollin,  and 
M.  A.  Cusanovich,  “Nonlinear  Optical  Properties  of  Photoactive  Yellow  Protein,” 
CLEO  '94. 

2.  Ljmgnes,  H.  M.  Gibbs,  C.  F.  Li,  S.  B.  Devanathan,  T.  E.  Meyer,  G.  Tollin,  and  M. 
A  Cusanovich,  “One-Photon  and  Two-Photon  Pump-Probe  Spectroscopy  of 
Photoactive  Yellow  Protein,”  J.  Nonlinear  Optical  Physics  and  Materials,  to  be 
published. 

3.  Lyngnes,  H.  M.  Gibbs,  C.  F.  Li,  S.  B.  Devanathan,  T.  E.  Meyer,  G.  Tollin,  and  M. 
A.  Cusanovich,  “One-Photon  and  Two-Photon  Pump-Probe  Spectroscopy  of 
Photoactive  Yellow  Protein,”  Final  Report  on  Rome  Laboratories 
F3060295C0100,  May  1996. 
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PIERRE  MEYSTRE’S  PUBLICATIONS 

1.  T.  Zaugg,  P.  Meystre,  G.  Lenz  and  M.  Wilkens,  “Theory  of  Adiabatic  Cooling  in 
Cavities,”  Phys.  Rev.  A  49,  301 1  (1994) . 

2.  P.  Meystre  and  M.  Wilkens,  “Spontaneous  Emission  by  Moving  Atoms,  in  Cavity 
Quantum  Electrodynamics,”  ed.  by  P.  R.  Berman  (Academic  Press,  New  York 
1994),  p.  301. 

3.  R.  Jin,  D.  Boggavarapu,  M.  Sargent  III,  P.  Meystre,  H.  M.  Gibbs  and  G.  Khitrova, 
“Photon-Number  Correlations  Near  the  Threshold  of  a  Microcavity  Laser,”  Phys. 
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Optics  VI,  ed.  by  J.  D.  Harvey  and  D.  F.  Walls,  (Springer,  Berlin  1994),  p.  68. 

5.  E.  Goldstein  and  P.  Meystre,  “Spontaneous  Emission  in  Optical  Cavities  —  a 
Tutorial  Review,”  in  Spontneous  Emission  and  Laser  Oscillation  in 
Microcavities,  ed.  by  H.  Yokoyama  and  K.  Ujihara  (CRC  press  1995),  p.  1. 
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8.  E.  Goldstein,  P.  Pax,  K.  J,  Schemthanner,  B.  Taylor  and  P.  Meystre,  “Influence  of 
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Trapping,”  Appl.  Phys.  B,  in  press. 
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NONLINEAR  OPTICS  AND  ELECTRON  OPTICS 
Galina  Khitrova 

Statement  of  the  Problem  Studied 

The  goal  was  to  grow  narrow-linewidth  quantum  wells  inside  a  semiconductor 
microcavity  by  MBE  and  study  the  linear,  nonlinear,  and  photoluminescence  properties 
of  this  strongly  coupled  system  of  exciton  and  cavity  mode.  Cavity  QED  effects  in 
semiconductor  microcavities  may  someday  lead  to  improved  lasers  and  light-emitting 
diodes. 

Abstract  of  the  Most  Important  Results 

InGaAs/GaAs  quantum  wells  with  record  narrow  linewidths  less  than  1  meV  have 
been  grown  inside  a  semiconductor  microcavity  consisting  of  GaAs/AlAs  distributed 
Bragg  reflectors  and  a  GaAs  spacer.  Vacuum  Rabi  splitting  or  normal  mode  coupling 
results  because  of  the  strong  coupling  between  these  two  oscillators,  namely  the  exciton 
and  the  cavity  mode.  The  quantum-well  narrow  linewidth  and  the  high-finesse  cavity 
result  in  record  values  of  the  splitting  to  linewidth  ratio.  This  permitted  the  observation 
of  a  new  nonlinear  effect:  with  increased  excitation  the  transmission  of  the  coupled 
system  goes  to  zero  with  negligible  reduction  in  the  splitting.  Measuremnts  on  quantum 
wells  outside  the  cavity  reveal  that  at  densities  below  the  onset  of  state  filling  and 
bandgap  renormalization  the  exciton  broadens  with  little  change  in  oscillator  strength. 
By  correlating  the  nonlinear  transmission  with  the  photoluminescence,  it  was  established 
that  the  density  is  too  high  for  excitons  to  behave  as  Bosons  when  the  upper  polariton 
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emission  increases  rapidly  compared  to  the  lower.  This  impressive  crossover,  called  by 
some  Boser  action  and  attributed  to  condensation  into  the  upper  polariton,  results  from 
the  density  dependence  of  both  the  emitters  and  the  photon  density  of  states. 

Summary  of  the  Principal  Findings 

Nonlinear  Behavior  of  Vacuum-Field  Rabi  Splitting  in  Semiconductor  Microcavities 
We  have  grown  InGaAs/GaAs  quantum  wells  with  very  narrow  linewidths  (1 
meV  =  0.6  nm).  By  growing  one  of  these  in  the  center  of  a  one-wavelength  GaAs  spacer 
or  two  in  the  two  central  anti-nodes  of  a  31/2  microcavity,  we  have  observed  record 
splitting-to-lineAvidth  ratios  (>10)  of  so  called  vacuum  Rabi  splitting  (VRS)  at  4K.  This 
phenomenon  is  also  known  as  normal  mode  coupling  because  it  exhibits  the  characteristic 
anti-crossing  curve  of  two  coupled  oscillators  (here  exciton  and  cavity)  as  their  relative 
detuning  is  scanned  through  zero.  We  have  increased  the  splitting,  seen  in  reflectivity  or 
transmission,  by  placing  the  microcavity  in  a  magnetic  field  up  to  12T,  showing  that  the 
splitting  depends  linearly  on  the  electric  dipole  moment  (published  in  Physical  Review 
B).  We  have  also  seen  record  depth  of  modulation  vacuum  field  Rabi  oscillations 
because  of  the  high  quality  of  our  samples.  We  have  studied  the  nonlinear  transmission 
and  reflection  of  our  VRS  microcavities  both  by  cw  pump/probe  and  by  femtosecond 
reflection  and  upconversion.  Because  of  our  very  narrow  lines  we  see  a  new  nonlinear 
behavior:  the  transmission  drops  to  zero  with  almost  no  reduction  in  VRS',  this  is 
explained  by  nonlinear  dispersion  theory  using  the  measured  exciton  absorption  with 
increased  carrier  density  in  a  transfer  matrix  formalism.  The  basic  physics  is  exciton 
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broadening  before  loss  of  oscillator  strength',  as  the  exciton  broadens  the  absorption 
increases  at  the  detuned  wavelengths  of  the  VRS  transmission  peaks  thus  decreasing  the 
Fabry-Perot  transmission  (a  little  bit  of  absorption  in  a  high-finesse  cavity  destroys  its 
transmission).  Since  the  area  under  the  exciton  absorption  (proportional  to  oscillator 
strength)  changes  very  little  as  the  broadening  sets  in,  there  is  almost  no  reduction  in 
splitting  while  the  transmission  drops  to  zero.  This  curious  nonlinear  behavior  could  not 
be  seen  by  other  groups  whose  VRS  samples  were  so  inhomogeneously  broadened  that 
exciton  broadening  had  little  effect.  The  cw  nonlinear  measurements  and  agreement  with 
a  many-body  theory  (with  Stephan  Koch  &  Frank  Jahnke,  Marburg,  Germany)  have  been 
published  in  Physical  Review  Letters.  Next  we  made  extensive  measurements  on  the 
photoliuninescence  spectra  of  VRS  samples  as  a  function  of  excitation  density  and 
correlated  them  with  the  nonlinear  transmission.  From  such  measurements  and  many- 
body  calculations  of  the  Bosonic  commutator,  we  established  that  the  rapid  increase  in 
the  upper  polariton  emission  relative  to  the  lower  polariton  is  not  a  result  of  condensation 
or  final  state  stimulation  as  proposed  by  others.  Instead  it  is  simply  the  emission  of  an 
excitation-dependent  distribution  of  emitters  within  a  microcavity  whose  photon  density 
of  states  is  also  excitation  dependent.  (Submitted  to  Phys.  Rev.  Lett.)  The  measurements 
described  so  far  were  performed  at  4K,  but  recently  we  have  achieved  record  room- 
temperature  normal-mode  coupling  using  InGaAs/GaAs  quantum  wells  in  a  microcavity 
with  aluminum-oxide/GaAs  mirrors  (published  in  Appl.  Phys.  Lett.). 
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PROTEIN 

Hyatt  M.  Gibbs  and  Michael  A.  Cusaonovich 

Statement  of  the  Problem  Studied 

Bacteriorhodopsin  (BR)  has  been  studied  and  promoted  as  a  medium  for  optical 
signal  processing.  Photoactive  yellow  protein  (PYP),  first  isolated  at  the  University  of 
Arizona,  is  a  smaller  and  better  characterized  molecule  with  linear  optical  properties 
similar  to  BR.  The  goal  of  this  research  project  was  to  measure  the  nonlinear  optical 
properties  of  PYP  and  to  evaluate  its  potential  for  optical  signal  processing  and  storage 
relative  to  PYP. 

Abstract  of  The  Most  Important  Results 

The  one-photon  nonlinear  absorption  of  PYP  was  measured  by  pump/probe 
spectroscopy.  Kramers-Kronig  transformations  on  those  data  yielded  the  corresponding 
nonlinear  refi'active  index  as  a  function  of  wavelength.  The  peak-wavelength  change  in 
index  per  excited  molecule  for  PYP  is  almost  the  same  as  that  of  PYP,  suggesting  that 
they  have  similar  one-photon  responses. 

An  extended  effort  was  made  to  measure  the  two-photon  absorption  cross  section 
of  PYP  without  success.  An  upper  limit  of  3.5  X  10*^2  s  molecule" ^  photon"  1  was 
determined,  10,000  times  smaller  than  that  reported  for  BR.  Consequently,  PYP  is  no 
competitor  with  BR  as  a  3D  optical  data  storage  medium,  if  the  reported  BR  cross 
section  is  correct.  Efforts  to  reproduce  the  BR  results  on  available  samples  were 
unsuccessful;  others  have  also  failed  to  reproduce  the  BR  results. 
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During  the  course  of  this  project  substantial  progress  was  made  in  the  methods  of 

growing  and  purifying  PYP  and  forming  new  forms  of  PYP  by  mutagenesis. 

The  details  of  this  entire  research  effort  are  contained  in  an  article  to  be  published 

in  the  Journal  of  Nonlinear  Optical  Physics  and  Materials. 
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ATOM  OPTICS 
P.  Meystre 
Research  Findings 

The  last  few  years  have  been  very  exciting  ones  indeed  in  the  fields  of  Quantum 
Optics  and  Atomic  Physics.  Many  exciting  new  developments  have  taken  place,  ranging 
from  quantum  state  preparation  and  control,  with  potential  applications  in  quantum 
computing,  to  atom  optics,  atomic  cooling  and  the  recent  realization  of  Bose 
condensation  in  alkali  vapors  and  the  first  demonstration  of  a  primitive  atom  laser. 

Our  earlier  research  and  research  plans  found  us  well  positioned  to  take  advantage 
of,  and  contribute  to,  these  latter  developments.  In  particular,  our  work  on  nonlinear 
atom  optics  gained  considerably  in  relevance  in  view  of  the  fact  that  low  density,  even 
multicomponent  quantum  degenerate  atomic  vapors  are  now  available. 

It  is  interesting  to  recall  that  in  our  1994  JSOP  annual  report,  we  stated: 

“...One  major  goal  for  the  next  years  is  to  fully  assess  the  potentiality  of  nonlinear  atom 
optics,  in  particular  in  the  realization  of  “coherent  atomic  beam  generators”,  sometimes 
called  “lasers  for  atoms.”  [...]  This  is  a  high-risk  field  of  investigations,  but  with 
considerable  potential  pay-off.  We  feel  that  the  realization  of  coherent  atomic  beam 
generators,  while  still  a  long-term  goal,  should  be  of  considerable  interest  both  for  the 
DoD  and  for  industry.  More  generally,  the  laser  manipulation  of  atomic  trajectories 
appears  to  rapidly  reach  the  level  of  sophistication  that  will  lead  to  emerging  new 
technologies  such  as  atomic  lithography,  etc... .” 
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At  that  time,  Bose-Einstein  condensation  of  low  density  atomic  vapors  was  still 
only  a  hope,  and  “atom  lasers”  sounded  almost  like  science  fiction.  In  retrospect,  the 
progress  achieved  since  then  has  been  nothing  short  of  amazing,  and  the  picture  of  an 
experimentally  realized  “atom  laser”  beam  even  made  the  cover  of  Physics  Today. 

Viewed  in  this  larger  context,  our  own  progress  has  been  rather  modest,  but  still, 
we  have  made  what  we  believe  to  be  a  series  of  useful,  and  sometimes  trailblazing 
contributions  to  that  field.  A  brief  summary  of  highlights  includes: 

Atom  Lasers:  Early  on,  we  have  proposed  an  atom  laser  scheme  using  the  conibined 
effects  of  dipole-dipole  interactions  and  Bose  enhancement  in  a  quasi-one-dimensional 
cavity  to  achieve  a  coherent  atomic  state  in  one  of  the  resonator  modes.  We  have  carried 
out  a  detailed  evaluation  of  the  dipole-dipole  selection  rules  between  center-of-mass 
modes  of  the  atoms  in  the  matter-waves  cavity,  and  studied  the  role  of  dimensionality  on 
these  cross-sections,  confirming  that  a  quasi  one-dimensional  cavity  appears  to  be  the 
geometry  of  choice. 

A  major  breakthrough  occurred  when  we  discovered  an  elegant  and  simple  way  to 
circumvent  the  effects  of  energy-momentum  mismatch  characteristic  of  most  light-matter 
wave  interactions:  While  the  dipole-dipole  matrix  elements  provide  rather  strict  center-of- 
mass  momentum  selection  rules,  these  are  largely  eliminated  by  the  energy  mismatch 
resulting  firom  the  quadratic  dispersion  relation  of  matter  waves.  We  found  that  this 
problem  can  be  completely  eliminated  by  modulating  the  light  fields  that  form  the  matter- 
waves  resonator  at  the  fi’equency  corresponding  to  this  energy  defect. 
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Following  this  work,  we  have  developed  a  master  equation  describing  the 
dynamics  of  the  atom  laser,  so  as  to  gain  an  understanding  of  the  dynamics  of  statistical 
properties  of  the  system,  its  linewidth,  etc.  We  have  determined  the  below-  and  above 
threshold  behaviors  and  identified  the  role  of  collisions  in  achieving  “laser  action.” 

Future  work  will  involve  the  study  of  higher-order  correlation  functions  of  the 
matter  field,  the  linewidth  of  the  atom  laser,  etc. 

Quantum  Transport  on  Optical  Lattices:  Another  aspect  of  nonlinear  atom  optics  that 
we  have  began  to  investigate  is  the  study  of  quantum  transport  in  optical  lattices:  The 
cooling  and  trapping  of  atoms  in  one-,  two-,  and  three-dimensional  optical  lattices  opens 
up  the  way  to  a  number  of  fascinating  investigations.  In  particular,  the  availability  of 
such  systems  provides  a  imique  probe  of  quantum  transport,  covering  situations  from  the 
extreme  case  where  deep  lattice  wells  tend  to  pin  the  atoms  to  sites,  so  that  their 
dynamics  becomes  analogous  to  that  of  a  spin  lattice,  to  the  full-fledged  quantum 
transport  of  atoms  between  wells.  Our  research  has  progressed  in  two  main  directions: 

First,  we  have  been  able  to  give  a  complete  static  description  of  the  near-resonant 
dipole-dipole  potential  in  3-dimensional  optical  lattices.  We  have  also  studied  this 
interaction  using  a  second-quantized  formalism  and  a  Wannier  fimctions  basis  to  evaluate 
the  effects  of  bound-boimd,  bound-fi-ee  and  fi-ee-fi-ee  transitions  on  lattices.  This  work 
will  be  useful  in  the  future  in  developing  a  full  manybody  theory  of  transport  on  lattices. 

Second,  we  have  carried  out  a  detailed  numerical  analysis  of  atomic  transport  in 
one-dimensional  lattices,  both  semiclassically  and  quantum  mechanically,  but  neglecting 


38 


so  far  the  effects  of  atom-atom  interactions.  Several  regimes  of  transport  have  been 
identified,  including  anomalous  diffusion  and  Levy  flights  for  intermediate  lattice  depths. 
While  there  are  quantitative  differences  between  the  quantum-mechanical  and 
semiclassical  predictions,  the  qualitative  results  of  these  two  approaches  are  rather 
similar.  This  is  an  encouraging  result  in  view  of  our  goal  to  study  transport  in  higher¬ 
dimensional  lattices,  where  full  quantum  calculations  are  beyond  the  capabilities  of 
present-day  computers. 

Nonlinear  Atom  Optics  and  Bose  Condensation:  In  addition  to  this  work,  we  have 
carried  out  research  more  directly  related  to  Bose  condensation  per  say.  For  example,  we 
have  studied  the  possibility  of  phase-conjugating  de  Broglie  waves  at  a  Bose  condensate. 
Both  our  theory  of  nonlinear  atom  optics,  as  well  as  the  Gross-Pitaevskii  description  of 
Bose  condensates  in  terms  of  a  nonlinear  Schrodinger  equation,  suggest  that  this  should 
be  the  case.  Phase  conjugation,  or  time  reversal,  of  matter  waves  would  be  of 
considerable  interest  both  in  the  context  of  fundamental  studies  of  irreversibility,  and  also 
as  possible  diagnostic  tools  for  Bose  condensates. 

In  addition,  we  have  begim  to  extend  the  Hartree-Fock  approach  to  nonlinear 
atom  optics  into  a  full  Hartree-Fock-Bogoliubov  theory.  In  contrast  to  the  case  of 
“traditional”  Bose  condensation,  which  involves  scalar  Schrodinger  fields  only,  general 
nonlinear  atom  optics  involves  vector  (multicomponent)  fields,  and  hence  leads  to  a  much 
richer  family  of  quasi-particle  excitations.  We  have  studied  the  dispersion  relations  of 
these  quasiparticles  in  some  simple  cases,  and  find  regimes  of  instability  where  the 
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condensate  solution  is  expected  to  be  subjected  to  modulational  instabilities.  Shortly  after 
our  first  results  were  sent  for  publication,  the  JILA  group  announced  the  experimental 
realization  of  a  two-component  condensate,  to  which  our  model  can  readily  be  applied. 
With  this  experimental  breakthrough,  the  study  of  multicomponent  matter  wave  fields  has 
become  reality  and  promises  to  become  very  important  in  the  near  future. 
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THEORY  OF  ULTRAFAST  NONLINEAR  OPTICAL  RESPONSE  OF 
UNIAXIALLY  STRAINED  SEMICONDUCTOR  QUANTUM  WELLS 

Rudolf  H.  Binder 

Abstract 

The  focus  of  this  project  was  twofold:  first,  linear  and  nonlinear  optical  properties 
of  uniaxially  strained  semiconductor  quantum  wells  have  been  investigated  theoretically 
and,  secondly,  a  theoretical  approach  to  the  vectorial  mode  structure  of  vertical-cavity 
surface-emitting  lasers  (VCSELs)  has  been  developed.  As  for  the  first  project,  a 
theoretical  description  for  the  ultrafast  (i.e.,  femtosecond  to  picosecond)  optical  response 
of  uniaxially  strained  GaAs  quantum  wells  has  been  developed.  This  theory  takes  into 
account  linear  and  coherent  nonlinear  exciton  effects  on  a  parameter-free  microscopic 
level.  Incoherent  effects  such  as  charge-carrier  scattering  were  considered  only  in  a 
preliminary  fashion  for  the  anisotropic  aspects  of  bulk  GaAs.  The  theoretical  results 
revealed  interesting  nonlinear  dynamics  of  the  interaction-induced  modification  of 
polarization  ellipticity  in  strained  quantum  wells.  As  for  the  vectorial  mode  structure  of 
VCSELs,  a  generalized  transfer  matrix  method  has  been  developed,  which  allows  for  a 
three-dimensional  vectorial  solution  of  Maxwell's  equations  for  micro-cavity  lasers  with 
cylindrical  symmetry. 
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Final  Report 

In  this  project,  we  have  investigated  theoretically  two  issues  from  the  general  area 
of  optical  communication  systems:  (i)  the  ultrafast  nonlinear  optical  response  of 
uniaxially  strained  semiconductor  quantum  wells  and  (ii)  the  vectorial  mode  structure  of 
vertical-cavity  surface-emitting  lasers  (VCSELs).  The  relation  between  these  two  issues 
lies  in  the  fact  that  VCSELs  do  not  necessarily  contain  perfect  cubic  quantum  wells  as 
active  material.  There  is  evidence  that  in  reality  these  quantum  wells  are  slightly 
strained,  a  generally  small  effect  which  can  have  a  large  effect  on  the  polarization 
properties  of  the  output  light  field  of  VCSELs.  A  solid  theoretical  foundation  was  laid 
for  both  issues  during  this  project,  and  we  anticipate  that  fiiture  projects  will  allow  us  to 
combine  both  theories  and  eventually  arrive  at  a  comprehensive  theoretical  model  of  the 
polarization  characteristics  of  micro-cavity  lasers  and  LEDs.  In  addition  to  these  long¬ 
term  goals,  the  investigation  of  the  nonlinear  optical  response  of  uniaxially  strained 
quantum  wells  complement  the  experimental  development  of  high-contrast  light 
modulators  at  the  Army  Research  Laboratory  (formerly  at  Fort  Monmouth).  Also  in  the 
context  of  microscopic  processes  in  anisotropic  structures  (ii),  we  have  investigated 
charge-carrier  scattering  and  quasi-thermalization  in  bulk  GaAs. 

All  of  these  subprojects  have  been  supported  through  the  JSOP  grant  as  well  as 
other  sources.  For  the  project  on  strained  quantum  wells,  JSOP  was  the  dominant  source 


of  funding. 
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The  first  issue  (i)  consists  of  several  subtopics,  out  of  which  the  coherent 
femtosecond  nonlinear  response  was  the  main  focus  of  this  project  and  required  the 
largest  amount  of  software  development. 

It  is  based  on  the  derivation  and  numerical  investigation  of  the  extended 
semiconductor  Bloch  equations  including  the  Luttinger  hamiltonian  to  account  for  strain 
effects  and  the  Coulomb  potential  to  account  for  linear  and  nonlinear  excitonic  effects. 
The  second  sub-topic  investigated  within  this  project  was  the  momentum-orientation 
relaxation  of  photo-excited  charge  carriers  by  means  of  carrier-carrier  scattering  in 
conventional  cubic  semiconductors,  such  as  ideal  bulk  GaAs.  This  laid  the  foundation  for 
future  investigations  of  carrier-carrier-scattering-induced  relaxation  and  quasi- 
thermalization  processes  in  anisotropically  strained  quantum  wells,  which  may  be 
combined  in  a  future  research  project  with  the  above-mentioned  theory  of  coherent 
nonlinear  optical  response  of  these  structures. 

All  of  the  above  listed  subprojects  have  lead  to  publications  in  peer-reviewed 
journals  [J1-J5]  and  were  presented  at  international  conferences  [C1-C8].  Ref.  [Cl]  was 
almost  finished  when  this  project  started.  It  is  a  comprehensive  review  article  on  the 
general  theory  of  the  semiconductor  Bloch  equations.  Ref.  [J2]  may  be  regarded  as  an 
additional  publication  in  the  field  of  femtosecond  optical  nonlinearities.  The  specific 
publication  regarding  the  semiconductor  Bloch  equations  for  anisotropic  systems  is  [J5] 
and  the  related  incoherent  scattering  theory  in  anisotropic  systems  is  [J3].  The  theory  of 
vectorial  electromagnetic  modes  of  VCSELs  is  presented  in  [J4]. 
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We  begin  with  the  description  of  the  theory  of  coherent  nonlinear  optical  response 
of  uniaxially  strained  GaAs  quantum  wells  [1,  2, 3]. 

The  theoretical  basis  is  the  Hartree-Fock  theory  of  optically  excited  thin 
uniaxially  strained  semiconductor  quantum  wells  which  involves  the  following  equations- 

of-motion  for  the  time-dependent  optical  polarization  functions  as  well  as  the 

electron  and  hole  distribution/coherence  functions  /„,(^)and  ,  respectively.  The 

vector  k  is  the  two-dimensional  in-plane  wave  vector,  the  electron  quantum  numbers 
s,s'  denote  the  spin-degenerate  conduction  bands  with  s  =  ±1/2,  and  the  hole  quantum 
numbers  j,j'  denote  the  two  degenerate  heavy-hole  (/  =  ±3/2)  and  light-hole  (/  =  ±1/2) 
bands.  The  equations-of-motion  (generalized  multiband  semiconductor  Bloch  equations 
[4]  are  ( ft  =1); 

%■/(*)! 

(1) 
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+  .  (2) 

s” 

and 

(^)  =  E  {f^M-  [^)f FT  (^)  -  fr  [^)^rr  (^) 

UL  j*» 

+ S  (n;-  (-  *)/^,(-  *)-«,(-  *)} 

5 

+S(S,.(‘)//5'(*)-4'(*)2:/,(*)}  ■  (3) 

r 

Here,  the  electron  energies  are  parabolic,  >  where  £g  is  the  effective 

band-gap.  The  Luttinger  hamiltonian  for  the  valence  bands  in  the  basis  { jj  =  (3/2,  -1/2, 
1/2,  -3/2)  with  basis  functions  defined  in  Ref  [5]  reads  (see  also  [1]) 


X  0  o' 

c  0  0 

0  0  c’ 

to  0  c 


(4) 


where  the  hh-hamiltonian  is  given 

by^AA  /2/MoX/', +r2)(^x +^j)-(l/2)A£;(l-2;',),  the  light-hole  (Ih) 
hamiltonian  is  /2m^Yi  +*j)  + A„..  +(l/2)A£;(l-2r.)and  the 
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hh-lh  coupling  is c  =  -(ft 2  / ‘^f’^o)^[r2 -K)~ ! 2)a£', .  The 
Luttinger  parameters  are,  y,,  yj,  yj,  and  the  strain  contributions 

-^w)*  with  =‘S',2 /(^ii +5,2),  contain  the  elastic 

compliances  5u  +  Sn ,  the  shear  deformation  potential  ft, ,  and  the  strain  tensor  elements 
e_„and  The  hydrostatic  stress  component,  which  merely  renormalizes  the  bandgap 

energy,  as  well  as  the  quantum  confinement  energy,  has  been  put  into  the  definition  of 
E'g  and,  therefore,  does  not  occur  explicitly.  The  hh-lh  splitting  is  denoted  byA^*_,*. 
The  energy  renormalizations  are  given  by  (with  a- sox  j) 

=  +  (5) 

and  the  renormalized  dipole  energy  is 


Cl^{k)  =  fi„-E  +  'Zv(q)P„{k*q). 

9 


(6) 


The  interaction  potential  v{q)  is  taken  to  be  the  unscreened  Coulomb  potential  in  two 
dimensions.  The  Jones  vector  of  the  negative-fi’equency  contribution  of  light  field, 
E  =  ,ey^ ,  can  be  written  as  [6] 
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E  =  (f  °  /  2%  +  (5;  /  2)  exp  (-  /0)e^  (7) 

=  [(^0  /  2)(cos^cos  £•  +  /  sin  0  sin  f)|e^ 

+  [(^’o  /  2)(sin dcos£-icos0 sin ,  (8) 

where  and  are  the  cartesian  unit  vectors,  s  is  the  ellipticity  angle  and  0  the  azimuth. 
Right-handed  circularly  polarized  light  is  defined  by  =  Sy  and  0  =  ^/2  and  linearly 
polarized  light  by  0  =  0.  The  dipole  matrix  elements  contain  the  information  of 

optical  selection  rules.  Using  circular  unit  vectors  one  can  write 

the  dipole  matrix  elements  [7]  as  f\n,in  -  ^nd 

M-V2,-V2  =  -1/2  =  where  ju  is  the  magnitude  of  the  microscopic  cartesian 

dipole  element.  Finally,  a  phenomenological  dephasing  contribution  with  dephasing  time 
T2  has  been  added  to  Eq.  1. 

From  a  numerical  point  of  view,  the  key  challenge  in  the  solution  of  the  above 
equations  was  the  evaluation  of  the  Coulomb  integrals,  which  are  of  the  form  of 
convolutions  of  the  Coulomb  potential  and  material  fimctions  such  as  the  polarization  and 

distribution  fimctions.  Since  ^  is  a  two-dimensional  vector,  the  integrals  are  two- 

dimensional  integrals  which  depend  on  the  two-dimensional  vector  k .  The  resulting 
numerical  complexity  had  been  the  reason  that  all  conventional  solutions  of  the 
semiconductor  Bloch  equations  had  been  restricted  to  isotropic  systems,  in  which  the 
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convolution  integrals  are  essentially  one-dimensional  and  depend  only  on  |^|.  The 
successful  completion  of  this  project  was,  therefore,  depending  on  the  question  whether 
optimized  algorithms  could  be  found  to  solve  the  anisotropic  semiconductor  Bloch 
equations.  It  turned  out  that  the  community  of  nuclear  physicists  had  successfully  dealt 
with  mathematically  similar  problems  by  means  of  Fast-Fourier-Transform  (FFT) 
algorithms.  The  basic  idea  is  to  use  FFT  to  transform  the  functions  to  be  convoluted  from 
momentum  into  real  space,  perform  the  product,  and  transform  the  product  back  via  FFT 
to  momentum  space.  This  yields  the  desired  convolution  integral.  The  application  of  this 
general  idea  to  the  specific  problem  of  the  multiband  semiconductor  Bloch  equations 
resulted  in  a  large  numerical  code  which  ran  on  Convex  as  well  as  Cray  C-90  computers. 
The  time-integration  is  performed  via  a  standard  4th  order  Runge-Kutta  method  and,  for 
the  calculation  of  optical  spectra,  the  subsequent  Fourier  transform  from  time  to 
frequency  space  is  performed  via  standard  integration  routines  (i.e.,  not  necessarily 
FFTs).  A  typical  time-integration  over  several  picoseconds  required  on  the  order  of  10 
minutes  to  1  hour  Cray  CPU  time. 

The  numerical  results  were  tested  thoroughly  to  yield  the  correct  behavior  for  the 
excitonic  dichroism  in  the  linear  optical  regime.  In  addition  to  comparisons  with 
analytically  known  limiting  cases  (i.e.,  linear  regime  without  Coulomb  interaction),  the 
comparison  with  experimental  results  obtained  at  ARL  confirmed  the  correct  use  of  the 
Luttinger  hamiltonian  and  the  exciton  contributions  to  the  equations-of-motion. 
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As  described  in  detail  in  Ref.  [J5],  very  interesting  results  can  be  obtained  in  the 
single-pulse  configuration  for  the  interaction-induced  polarization  rotation.  Here,  the 
repulsive  interaction  between  excitons  is  the  main  cause  of  the  nonlinear  response.  Pauli- 
blocking  and  saturation  effects  have  been  found  to  be  less  dramatic. 

As  for  the  issue  of  incoherent  scattering  processes,  the  corresponding  subproject 
was  integrated  into  a  collaboration  with  a  group  of  nuclear  physicists.  The  results  on 
momentum-orientation  relaxation  via  electron-electron  scattering  in  cubic  bulk  GaAs 
crystals  was  published  in  Ref  [J3].  The  theoretical  basis  is  the  incoherent  part  of  the 
generalized  semiconductor  Bloch  equations.  This  part  is  best  known  as  generalized 
quantum  Boltzmann  equations  or  Kadanoff-Baym  equations.  For  a  one-component 

electron  plasma,  these  equations  describe  the  time-evolution  of  the  two-time  one-particle 

> 

Green's  fimctions 

^0 

/'  > 

^0 


(9) 


54 


Within  the  screened-Hartree  Fock  approximation  and  the  additional  approximation  of 

> 

quasi-static  screening,  the  self  energy  Z^(^,r,r')is  given  by 

>  >  >  > 

k\q 

where  W{^,t)  =  ^ne^  +^^(0)  the  quasi-statically  screened  Coulomb  potential  in 

Gauss  units.  Here,  =  el!  4  is  the  effective  charge  squared  is  the  electron  charge  in 
vacuum  andf^is  the  background  dielectric  constant  of  the  semiconductor)  and  k  is  the 
screening  wavenumber.  As  for  the  initial  condition,  we  assume  that  at  the  initial  time 

/gall  correlations  vanish  and  that  the  equal- time  Green's  function  the 

initial  distribution  discussed  above.  Due  to  the  specific  retardation  described  by  Eqs.  9 
and  10,  the  initial  condition  at  a  single  point  in  the  t-t'  plane  is  sufficient  to  obtain  the 
solution  for  all  t,t'>tQ.  The  carrier  distribution  function  is  finally  obtained  from 

The  basic  development  of  numerical  algorithms  of  the  anisotropic  version  of  the 
Kadanoff-Baym  equations  was  not  part  of  this  project,  for  they  had  been  developed  by 
the  nuclear  physicist  taking  part  in  this  collaboration.  The  main  goal  of  this  project  was 
the  application  of  the  Kadanoff-Baym  equations  to  the  case  of  a  photo-excited  electron 
plasma  in  a  semiconductor.  The  hole  plasma  has  been  neglected  to  simplify  the  problem. 
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To  include  the  effects  based  on  the  Luttinger  hamiltonian  into  the  analysis,  the 
initial  condition  of  the  time  integration  of  the  Kadanoff-Baym  equations  were  chosen  in 
the  following  way.  We  followed  Zakharchenya  et  al.  [8]  and  used  an  angle  dependence 
of  the  two  electron  components  (one  due  to  hh-c  excitation  and  the  other  due  to  Ih-c 
excitation)  that  results  form  perturbation  theory  of  the  optical  excitation  process  with 

linearly  polarized  light.  We  denote  the  angle  between  a  given  momentum  state  k  and  the 
polarization  vector  of  the  light  as  0.  As  shown  in  [8]  the  momentum  dependent  squared 

dipole  matrix  elements  |<  c,^|r -.Elv.A  >|^  between  the  conduction  band  “c”  and  the 
valence  band  “v”,  where  v  stands  either  for  hh  or  Ih  denoting  heavy-hole  and  light-hole, 
respectively,  is  proportional  to  IT P2(cos0) ,  where  gives  the  contribution  due  to  the 
hh-c  transition  and  “+”  that  of  the  Ih-c  transition.  Pj  is  the  second  degree  Legendre 
polynomial.  Not  given  in  [8]  are  the  proportionality  factors  which  depend  only  on  the 
modulus  of  k  and  the  details  of  the  optical  excitation  such  as  the  amplitude,  duration, 

and  center-frequency  of  the  optical  excitation  field  E .  In  our  analysis,  we  proceeded 
quasi-phenomenologically  and  used  prefactors  which  correspond  to  an  optical  excitation 
process  in  which  the  distribution  function  reflects  the  spectral  shape  of  the  optical  pulse 
which  we  take  to  be  Gaussian.  In  the  context  of  a  third-order  approximation  of  the 
semiconductor  Bloch  equations  without  Coulomb  interaction,  this  can  be  achieved  in  an 
approximate  way,  if,  in  the  time-derivative  of  the  distribution  function,  the  interband 
polarization  function  is  always  taken  at  times  long  after  the  pulse.  This  procedure  yields 
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Gaussian  distribution  functions  centered  at  the  carrier  momentum  that  corresponds  to  the 
resonance  condition  of  the  excitation  field:  HcOq  =  +s°  with  a  =  Here,  (Vg  is 

the  center  frequency  of  the  light  pulse,  1 2m^  +  Eq  is  the  electron  energy  with 

,,  ,  V 

effective  mass  is  the  bandgap  energy,  =— — [y^  +2/3] is  the  heavy-hole 

ZJfto 

energy  with  Luttinger  parameters  ^2  ("*o  the  electron  mass  in  vacuum)  and 
Sr  = - (y.  +  ly. )  is  the  light-hole  energy.  Our  initial  distribution,  therefore, 

*  Intg  ^  ' 

contains  two  contributions  with  slightly  different  center-momenta  (due  to  the  mass 
difference  of  hh  and  Ih)  and  different  angle  characteristics  because  the  hh  has  an  angular 
momentum  quantum  number  of  7  =  ±3  /  2  (with  respect  to  the  quantization  axis  given  by 

k  ),  whereas  the  Ih  has  7  =  ±1  /  2 : 


/(^,0)  =  A^^{k)^\\ -  cos^  0] 
+  ^/a(^)7[i  +  3cos^  0] 


(11) 


with 


4  (^)  =  Ag  exp(-  [hcOg  -  si  -  el  ^  /  (ly^fj , 


(12) 


/  \2  /r 

where  Ag=yiLigSgj  'TT*  A)  Is  the  angle-independent  part  of  the  interband  dipole 


moment,  ^g  is  the  peak  light-field  amplitude  and  the  width  of  the  distribution  is  related  to 
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the  pulse  duration  t  (FWHM  of  intensity)  through;'^  =^^21n2/x^.  This  approach,  in 
which  the  light-field  parameters  enter  only  into  the  shape  of  the  initial  distribution,  is 
strictly  justified  only  if  the  relaxation  processes  are  much  slower  than  the  light-pulse 
duration  x. 

The  numerical  solutions  of  the  initial  value  problem  were  obtained  with 
algorithms  similar  to  the  ones  discussed  above  in  the  context  of  the  coherent  anisotropic 
semiconductor  Bloch  equations.  The  algorithms  were  the  same  as  those  used  in  Ref.  [9]. 
The  maximum  k-value  (in  each  Cartesian  direction)  was  chosen  to  be  between  6.3  and  8.7 

in  units  of  Og'  (where  Ug  is  the  exciton  Bohr  radius  of  the  systems,  which  is  here 

o 

Og  =  132  A)  and  the  number  of  k-points  in  each  direction  was  between  43  and  59. 
Within  this  parameter  range,  the  results  were  foimd  to  be  insensitive  to  the  parameters. 

The  solutions  revealed  the  unexpected  result  that  the  anisotropy,  measured  by  the 
quadruple  moments  of  the  distribution  function  in  momentum  space,  decays  slower  than 

the  zone-center  occupation /(O,/),  which  is  a  measure  for  the  monotonicity  of  the 
distribution  function  in  momentum  space,  increases.  Further  discussion  can  be  found  in 
Ref.  [J3]. 

Finally,  we  turn  to  the  topic  which  was  designed  to  lay  the  foundation  for  a  future 
application  of  the  microscopic  theories  discussed  above  to  the  technology-oriented  field 
of  micro-cavity  lasers  and  VCSELs.  This  work  was  partly  motivated  by  the  fact  that  the 
understanding  of  the  polarization  characteristics  of  the  VCSEL  output  requires  a  detailed 
knowledge  of  the  vectorial  eigenmodes  of  the  empty  cavity.  The  understanding  of  the 
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vectorial  mode  structure  provides  the  proper  basis  in  which  one  can  expand  the  light 
field  when  modeling  running  VCSELs.  This,  in  turn,  should  enable  the  understanding  of 
the  experimentally  observed  coexistence  and  switching  between  transverse  modes  [10- 
18]. 

One  goal  of  this  project  was  to  provide  a  general  as  well  as  an  approximate 
analytical  method  to  analyze  the  modal  structure  of  VCSELs  taking  into  account  their 
polarization  properties.  For  clarity,  we  focused  on  the  idealized  system  of  cold-cavity 
modes,  i.  e.,  we  assumed  all  refiractive  indices  to  be  real,  and,  therefore,  the  absence  of 
gain  and  absorption  fi’om  the  system.  In  the  cold-cavity  configuration,  the  only 
information  about  the  electromagnetic  modes  of  VCSEL  structures  can  be  obtained  fi’om 
their  reflectivity  spectra,  which  was  also  calculated. 

The  general  approach  was  to  apply  the  knowledge  of  electromagnetic  hybrid 
modes  in  optical  fibers  and,  based  on  these  modes,  develop  a  generalized  vectorial 
transform  matrix  method  to  allow  for  a  three-dimensional  solution  of  Maxwell's  equation 
in  cylindrical  symmetry.  The  starting  point,  therefore,  is 

V  X  £  =  and 


c 


(13) 
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where  c- 1/ and  //oare  speed  of  light,  electric  and  magnetic  field 
permeabilities  of  a  fi-ee  space,  respectively,  n,  is  the  refractive  index  of  the 
corresponding  dielectric  layer,  f  =  0,l,..A/^  +  l.  In  cylindrical  coordinates,  the  field 

components  are  E  =  {e^  E^E.^  and  H  =  We  assume  that  the  solution  of 

Eqs.  (13)  is  of  the  general  form  £(r)0(z,^)and  H{r)2^z,^) ,  where 

0(z,^)  =  ,  with  being  the  z  -  component  of  the  wave  vector,  and  m  being  an 

integer  describing  the  mode  properties  in  azimuthal  direction.  The  general  solutions  can 
be  found  in  textbooks  on  fiber  optics,  and  they  are  also  listed  in  Ref.  [J4].  The  conditions 
of  continuity  at  the  VCSEL  radius  r=R  yields  the  well-known  dispersion  relation  of  the 
hybird  modes 


"1  j'.  ,  ^  k] 

^  ^  Jin  1  ^ 

r  1  i] 

J„  Kj 

W  Kj 

~  R^kl^ 

where  (/K)  and  (SR)  are  Bessel  functions  evaluated  at  the  boundary, 

and  the  transverse  propagation  constants  are  given  by 


(15) 


inside  and 


S  =  ^k‘-k’. 


(16) 
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outside  the  active  layer  of  the  VCSEL.  The  active  layer  has  the  layer  number 
=  IN ^  + 1  and  we  denote  =  P  and  • 

The  transverse  components  of  the  electric  and  magnetic  field  in  the  ^  -  th  layer 
may  be  expressed  as  a  superposition  in  a  complete  set  of  orthogonal  modes  of  the 
corresponding  cylindrical  waveguide  as 

s',  =  lLs’^  and  =  (17) 

V  y 

with 

e'^  =  (a,'V»''  -  (18) 

it'^  (19) 

where  the  subscript ‘T  ”  stands  for  transverse  field  component.  Here,  we  have  explicitly 
written  the  forward  and  backward  waves  and  the  vectors  ^^(T,^)and  ^^(T,^)are 

defined  by  comparing  Eqs.  (17)  with  the  full  definition  of  the  transverse  field  components 
given  in  Table  I  and  n  of  Ref.  [J4].  Using  the  appropriate  boundary  conditions  at  the 
interface  between  two  layers  (say  i  and  ^  + 1 ),  we  obtain  the  final  vectorial  transform 
matrix  formulation,  which  relates  the  amplitudes  of  the  forward  and  backward  running 
waves  of  each  hybrid  mode  in  adjacent  layers. 
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^e+i 


=  Mi 


'5' 


^  Mr  m; 

with  M(  = 

Mr  Mr 


(20) 


where  M<  is  the  transform  matrix  at  the  z  =  z, . 

In  order  to  simplify  the  notation  in  Eq.  (20),  we  introduced  vectors  (A*  and 

A^^,)  and  matrices  ( M^ )  with  elements  given  by 


A^=({Ar}y,A^.,=i{A^:,}y, 


(21) 


-p- 

rr.±r: 


(22) 


a+\ 


->(ri^ri*i)2e 


YM+Yt 


(23) 


where  the  superscript ‘T”  indicates  transposed  vector  and 


Q 


tM 

VM 


i  Efr  (  2  _  2)  \A,dA,El(v,(p)-El*'(T,(p) 

2V>“o  ”  '  \^dA^E‘;'(t,(p)xH‘;'(r,<p)  -f 


(24) 
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with  Af  denoting  the  transverse  cross-section  at  the  z  =  Zf  plane. 

Finally,  multiplication  of  the  transform  matrices  Mi  for  all  interfaces  leads  to  the 
total  transfer  matrix  of  the  VCSEL  structure 

_  N  _ 

M  =  YlMt.  (25) 

;=o 


In  order  to  determine  the  eigenmodes  of  the  VCSEL  cavity,  we  note  that  in  the 
running  laser  only  the  outcoming  amplitudes  are  nonzero,  that  is  =  0  and  A^  =  0. 

Therefore,  the  condition  for  the  cavity  eigenmode  has  the  form:  M~Al  =  0  or, 

equivalently,  det  M~~  =0,  where  Af''is  a  submatrix  of  the  total  matrix  M  (see  Eq. 
20).  The  solution  of  this  condition  yields  resonant  wavelengths,  threshold  gain  and 
eigenmode  profiles  for  the  VCSEL  cavity. 

There  are  great  advantages  to  such  an  almost  analytical  solution  of  Maxwell's 
equations,  especially  for  the  future  project  of  developing  a  comprehensive  microscopic 
model  of  the  laser  operation  of  VCSELs.  Since  the  microscopic  models  for  the  nonlinear 
optical  behavior  and  optical  gain  characteristics  of  the  active  layers  require  extensive 
numerical  resources,  the  semi-analytical  solution  of  the  corresponding  electromagnetic 
problem  is  needed  to  keep  the  problem  numerically  feasible.  But  the  use  and  applications 
of  the  vectorial  transform  matrix  approach  does  not  require  the  targeted  full  microscopic 
model.  Much  information  about  the  dependence  of  laser  threshold  and  resonant 
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wavelengths  can  already  be  obtained  from  the  electromagnetic  model  alone,  as  discussed 
for  example  in  Ref.  [J4]. 
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HIGH  POWER  CIRCULAR-GRATING  SURFACE-EMITTING  LASERS 

Mahmoud  Fallahi 

Objectives: 

Investigate  issues  to  the  high  power  operation  of  circular-grating  surface-emitting 
DBR  (CG-SE-DBR)  lasers. 

The  effect  of  grating  positioning  and  aspect  ratio  on  the  lasing  performance  (threshold 
and  output  power). 

Fabrication  of  high  power  CG-SE-DBR  lasers  on  InGaAs/GaAs  multi-quantum-well 
structure. 

Status  of  Effort: 

Circular-grating  surface-emitting  DBR  lasers  are  very  attractive  for  high  power 
application.  Because  of  their  two  dimensional  nature,  the  position  control  and  uniformity 
of  the  circular  gratings  are  the  key  factors  in  low  threshold  high  power  operation  of  the 
lasers.  The  inclusion  of  a  thin  GaAs  etch-stop  layer  in  the  top  AlGaAs  cladding  layer  for 
position  control  resulted  in  the  lowest  threshold  ciurent  for  CG-SE-\DBR  lasers  ever 
reported.  A  threshold  current  less  than  20  mA  and  a  pulsed  output  power  in  excess  of 
170  mW  is  successfully  demonstrated.  This  also  resulted  in  a  reproducible  operation  of 
the  lasers  across  the  wafer.  The  laser  operated  under  CW  pumping  with  the  same 
threshold  current  but  the  output  power  saturated  around  7.5  mW  caused  by  a  poor  heat 
sink. 

Technical  Summary: 

There  is  a  growing  interest  in  low  divergence,  stable  single-mode,  high-power 
surface-emitting  lasers  for  various  applications.  Among  the  potential  devices,  circular- 
grating  surface-emitting  DBR  lasers  are  attractive  candidates  to  provide  single  mode, 
low-divergence,  high-power  surface  emission.  Because  of  the  two  dimensional  nature  of 
these  lasers,  the  position  control  and  the  imiformity  of  circular  gratings  are  the  key 
factors.  The  addition  of  an  etch-stop  layer  in  a  strained  multiple  quantum  wells  (MQW) 
structure  is  an  efficient  solution  allowing  both  the  use  of  the  same  structure  for  the  active 
and  passive  sections  as  well  as  an  accurate  positioning  for  the  grating.  In  this  work  we 
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successfully  fabricated  low  threshold  current  and  high-power  operation  of  an  electrically 
pumped  CG-SE-DBR  laser  using  a  strained  triple-quantum-well  (TQW)  structure  with  an 
etch-stop  layer. 

Fig.  1  shows  a  schematic  sectional  view  of  the  device.  Similar  to  the  case  of 
linear  grating-coupled  DBR  lasers,  the  two  major  parameters  to  be  considered  in  the 
design  of  the  CG-DBR  lasers  are  the  power  coupling  efficiency  between  the  gain  region 
and  the  grating  region  (Cq)  and  the  coupling  coefficient  (k)  between  the  inward-  and 
outward  propagating  waves.  For  a  four  layer  structure  they  can  be  defined  as  follow  [5]: 


and 


Co  = 


\ZcZgdz 


(1) 


k  =  T, 


\.2J3 


Sclad  —  €  cover 


J 


wn 


Wg 

sin  nm  — . 
V  AJ 


(2) 


where  Za  and  Zg  are  the  field  distribution  for  the  active  and  grating  sections  in  the  z 
direction,  ko  is  the  wave  number  of  light  in  vacuum,  P  is  the  propagation  constant  of  the 
unperturbed  waveguide,  m  is  the  grating  order,  sdad  and  ^cover  are  the  permittivity  of 
the  cladding  layer  and  the  cover  layer  (air  in  our  case)  respectively,  wg  is  the  width  of  the 
grating  valley.  Fg  is  the  overlap  integral  given  by: 

u=ljzfdz  iwn'dz  (3) 

The  eigenvalue  equation  for  the  threshold  condition  can  be  approximated  as: 

CopopKxe  e  =  1  (4) 

Po  is  the  reflection  coefficient  at  r  =  0  and  is  defined  as  the  ratio  between  outward- 
propagating  to  the  inward-propagating  waves,  is  the  field  reflection  coefficient  at  r  = 
Rl.  a  represents  the  gain  (or  loss)  and  Q  is  the  total  phase  shift  given  by: 
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^  Tn{wg+  2i?l)+  l^ij^gain  —  ^grat^Rl  (^) 

For  a  CG-DBR  as  presented  in  fig.l,  we  have  po  =  1  while  pRi  depends  strongly  on  the 
coupling  coefficient  k  as  follow: 

sinh(?£) 

^  )'cosh{yL)-{ag-JS)siDh(/L) 

with 

r  =  +(ag-jsy  (7) 

where  5  is  the  deviation  from  the  Bragg  frequency,  L  is  the  grating  length,  n  is  the  order 
of  the  cylindrical  waves  and  ag  is  the  loss  in  the  grating  region. 

The  above  equations  show  that  in  order  to  reduce  the  threshold  current,  high 
coupling  efficiency  and  coupling  coefficient  are  needed.  The  use  of  the  same  layer 
structure  for  the  gain  section  and  grating  section  is  a  simple  method  of  obtaining  high 
coupling  and  eliminating  the  need  for  epitaxial  regrowth.  In  this  case,  the  value  of 
coupling  coefficient  and  power  coupling  efficiency  are  strongly  influenced  by  the 
position  of  the  grating  (and  cladding  thickness). 

Circular-grating  surface-emitting  DBR  lasers  similar  to  Fig.  1  are  designed  and 
fabricated.  The  layer  structure  is  an  InGaAs/GaAs  strained  TQW,  graded-index,  separate- 
confinement  heterostructure  (GRINSCH).  A  GaAs  etch-stop  layer  is  incorporated  in  the 
top  cladding  structme.  The  circular  pumped  region  of  the  laser  is  placed  in  the  center  of 
the  structure  while  the  passive  Bragg  section  surroimds  the  active  section.  Second  order 
circular  gratings  with  a  radial  extension  of  150  pm  were  defined  by  an  electron-beam 
system  with  polar  coordinate  writing  capability. 

The  structure  was  grown  by  molecular  beam  epitaxy  (MBE)  on  an  n'^'-GaAs 
substrate.  The  active  region  consists  of  three  wells  of  60  A  of  frixGai-xAs  (x=20  %) 
separated  by  300  A  GaAs  barriers.  This  active  region  was  sandwiched  by  graded  AlxGai- 
xAs  layers.  Top  and  bottom  cladding  consist  of  AlxGai-xAs  (x=67%)  layers.  A  150  A 
GaAs  etch-stop  layer  was  grown  in  the  top  cladding  layer,  0.2  pm  above  the  graded 
region.  A  60  pm  diameter  pumping  region  was  defined  by  reactive  ion  etching  (RIE)  of 
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the  top  cladding  layer  in  BCl3:He  down  to  about  0.2  pm  above  the  etch-stop  layer.  Then 
the  remaining  AlGaAs  was  selectively  etched  down  to  the  etch-stop  layer  in  HF  solution. 
Next  500  A  of  Si02  was  evaporated  on  the  structure. 


MQW 


Fig.  1  :  Schematic  view  of  a  CG-SE-DBR  laser 

Concentric  circular  gratings  were  defined  by  electron-beam  lithography  in  0.25 
pm  PMMA  resist  [6].  The  radial  extent  of  the  circular  gratings  surrounding  the  active 
region  is  about  150  pm.  The  grating  was  next  transferred  into  the  oxide  by  RIE  in 
CF4:02  .  The  gratings,  patterned  in  the  oxide,  were  then  transferred  into  AlGaAs  by  RIE 
in  BCl3:He  using  oxide  as  a  mask.  Since  the  oxide  is  imafifected  by  the  etch,  a  high  aspect 
ratio  and  an  optimized  coupling  coefficient  can  be  achieved.  Finally  an  n-contact  was 
evaporated  at  the  back  of  the  sample  and  the  contacts  were  annealed. 

The  measurements  were  done  at  room  temperature.  The  lasers  were  characterized 
under  pulsed  and  continuous  wave  (cw)  operation.  The  threshold  current  is  less  than  20 
mA  for  both  cases.  The  external  quantum  efficiency  was  as  high  as  14  %.  An  output 
power  in  excess  of  170  mW  without  saturation  under  pulsed  operation  was  obtained.  The 
cw  operation  showed  no  change  in  the  threshold  current  but  the  output  power  saturated  at 
about  7.5  mW  for  a  150  mA  current  most  probably  due  to  the  relatively  high  resistivity  of 
the  p-cladding  layer  and  poor  heat-sinking  The  lasing  wavelength  was  measured  to  be 
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about  980  run.  Fig,  2  shows  the  lasing  spectrum  of  a  CG-SE-DBR  laser  at  2  times  the 
threshold  current.  Single-mode  operation  at  low  injection  and  a  mode  hopping  at  high 
injection  were  observed. 


\/\^ength(niTj) 


Fig.  2 :  Lasing 


spectrum  of  a  CG-SE-DBR  laser  at  21* 


We  measured  a  farfield  divergence  angle  of  less  than  2  degrees  in  all  directions  with  a  quasi-circular 

cross  section.  At  high-power  operation,  some  of  the  lasers  showed  a  multi-lobe  far  field  pattern 
with  a  total  divergence  angle  as  high  as  4.S  degrees,  due  to  the  multi-mode  operation  of  those 
lasers.  These  results  indicate  the  great  potential  of  circular-grating  surface-emitting  DBR  lasers 
for  high  power  applications.  In  addition  since  no  cleaved-fact  is  required,  a  2-D  array  of  these 
lasers  can  be  fabricated. 
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ATOM  TRAPPING  IN  FAR-OFF-RESONANT  OPTICAL  LATTICES 

Poul  Jessen 

Introduction. 

The  AC  Stark  shift  (light  shift)  arising  in  laser  interference  patterns  can  be  used  to 
create  stable  periodic  potentials  for  neutral  atoms.  Under  appropriate  conditions  these 
"optical  lattices"  will  laser  cool  and  trap  atoms  in  individual  optical  potential  wells,  with 
center-of-mass  motion  in  the  quantum  regime.  Experiments  have  explored  a  large 
number  of  lattice  configurations  in  one,  two  and  three  dimensions,  and  have  used  a 
variety  of  techniques,  including  probe  absorption  and  fluorescence  spectroscopy,  phase 
conjugation  and  Bragg  scattering,  to  gain  detailed  insight  into  the  dynamics  of  cooling 
and  trapping  [1].  In  particular,  it  has  been  demonstrated  that  atoms  are  readily  localized 
deep  in  the  Lamb-Dicke  regime.  This  suggests  that  it  is  worthwhile  to  pursue  resolved- 
sideband  Raman  cooling  [2]  and  quantum  state  preparation  [3],  by  employing  techniques 
similar  to  those  recently  demonstrated  for  trapped  ions. 

In  an  optical  lattice  formed  by  near-resonance  light,  control  of  the  center-of-mass 
motion  is  limited  by  laser  cooling  and  heating  processes  that  occur  at  a  rate  determined 
by  photon  scattering.  These  dissipative  processes  can  be  avoided  if  the  lattice  is  formed 
by  intense  light  tuned  far  from  atomic  transition.  Such  far-off-resonance  lattices  have 
been  used  extensively  in  atom  optics  as  diffraction  gratings  and  lenses  [4],  and  recently  as 
model  systems  in  which  to  study  quantum  chaos  [5]  and  quantum  transport  [6].  The 
research  supported  by  this  grant  has  studied  the  use  of  far-off-resonance  optical  lattices  to 
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trap  neutral  atoms  in  the  Lamb-Dicke  and  quantum  regime,  and  in  an  environment 
essentially  free  of  dissipation.  A  key  goal  has  been  to  evaluate  the  feasibility  of  resolved- 
sideband  Raman  cooling  and  coherent  control  of  the  atomic  center-of-mass  motion. 
Summary  of  Research 

1.  Demonstration  of  a  far-off-resonance  optical  lattice  for  Cesium  atoms 

An  important  accomplishment,  and  a  prerequisite  for  much  of  the  research 
discussed  below,  was  our  demonstration  of  Cesium  atom  trapping  in  a  1 -dimensional 
(ID)  far-off-resonance  optical  lattice.  The  lattice  was  formed  using  a  0.5  W 
semiconductor  laser  tuned  anywhere  between  5  and  40  GHz  below  the  Cs 
6Sf2  (^  =  4)-^  6/^32  (j^'  =  5)  transition  at  A  =  852  nm . 

2.  Loading  scheme  for  deeply  boimd  states  of  a  far-off-resonance  optical  lattice. 

When  atoms  are  trapped  in  a  far-off-resonance  optical  lattice,  the  absence  of  built- 
in  laser  cooling  makes  it  difficult  to  obtain  vibrational  excitation  and  confinement 
comparable  to  the  near-resonance  case  [7].  This  in  turn  makes  sideband  cooling  schemes, 
as  well  as  related  means  of  coherent  control,  considerably  less  efficient.  A  key  part  of 
our  research  has  therefore  been  the  development  of  a  loading  scheme  to  overcome  this 
limitation  [8].  Our  technique  relies  on  adiabatic  transfer  of  Cesium  atoms,  initially 
cooled  and  trapped  in  a  near-resonance  lattice,  into  a  superimposed  far-off-resonance 
lattice.  If  truly  adiabatic  transfer  can  be  accomplished,  then  the  net  result  will  be  a  one- 
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Fig.  1.  (a)  Mean  vibrational  excitation  n  in  the  far-ofF-resonance  lattice,  as  a 
function  of  the  peak-peak  modulation  depth  Upp  of  the  lattice  potential.  Solid  (open) 
symbols  indicate  data  taken  for  lattice  detunings  of  A  =  -20  GHz  (A  =  -10  GHz). 
Circles  (squares)  indicate  vibrational  excitation  measured  at  r=  20  ms  (r=  0.1  ms) 
after  transfer  fi^om  the  near-resonance  lattice,  (b)  Lamb-Dicke  parameter  (kAz)^ 
corresponding  to  the  data  in  (a). 

transfer  of  atoms  between  bound  states  of  the  two  lattices.  Immediately  following 
transfer  we  do  in  fact  achieve  trapping  parameters  comparable  to  the  near-resonance  case, 
with  a  mean  vibrational  excitation  as  low  as  n  «  0.3  (consistent  with  no  vibrational 
heating  during  transfer),  and  a  typical  rms  position  spread  of  Az  «  A/20  (fig.  1).  Transfer 
efficiency  between  the  near-and  far-off-resonance  lattices  is  typically  in  the  range  90- 
95%.  A  small  bias  field  applied  while  atoms  are  cooled  and  trapped  in  the  near-resonant 
lattice  allows  us  to  skew  the  distribution  of  population  towards  one  of  the  magnetic 
sublevels  \m  =  ±4) ,  as  discussed  in  more  detail  below.  A  small  bias  magnetic  field, 
applied  during  transfer,  helps  maintain  the  atomic  orientation.  With  those  minor 
modifications  of  the  loading  scheme  we  can  preferentially  populate  a  specific  magnetic- 
and  vibrational  groimd  state  of  the  far-off-resonance  lattice. 
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3.  Trapping  dynamics  in  a  far-off-resonance  optical  lattice. 

Photon  scattering  constitutes  a  major  source  of  heating  for  atoms  trapped  in  an 
optical  lattice.  Fig.  2  shows  the  increase  in  mean  vibrational  excitation  «  ,  as  a  function 
of  trapping  time  r,  measured  for  lattices  of  different  depth  and  detuning  [8].  Initially  n 
increases  roughly  as  expected  from  photon  scattering;  fig.  2  shows  calculations  of  n  vs. 
time,  obtained  by  solving  rate  equations  for  the  vibrational  populations  for  an  atom  bound 
in  an  optical  potential  well  on  the  axis  of  the  lattice  beams.  Note  that  one  would  expect 
the  rate  of  heating  to  increase  with  photon  scattering  rate,  a  behavior  which  we  observe 
only  during  the  first  few  ms.  At  later  times  n  depends  only  on  the  lattice  depth.  Fig.  la 
illustrates  this  linear  scaling  of  n  with  lattice  depth,  at  r=20ms.  The  observed 
behavior  is  qualitatively  similar  to  the  steady  state  scaling  known  from  near-resonance 
lattices  [1],  and  one  might  speculate  that  a  cooling  mechanism  is  active  also  in  the  far-off- 


Fig.  2.  Mean  vibrational  excitation  n  in  the  far-off-resonance  lattice,  as  a  function  of  time  r 
elapsed  since  transfer  from  the  near-resonance  lattice.  Solid  circles  correspond  to  a  lattice 
detuning  A  =  -20  GHz  (-3831  F)  and  a  peak-peak  modulation  depth  Upp  =  of  the  potential. 
Open  symbols  correspond  to  A  =  10  GHz  (-1916  T)  and  Upp  =  105  £j;  (circles),  Upp  =  209  Er 
(squares).  Solid,  dashed  and  dot-dashed  lines  show  the  expected  heating  from  photon  scattering. 
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resonance  case.  Unfortunately  the  measurement  uncertainties  and  limited  interaction 
time  imposed  by  our  ID  geometry  does  not  permit  us  to  confirm  or  rule  out  whether  a 
steady  state  will  eventually  be  reached. 

It  is  possible  that  the  increase  in  n  is  limited  by  the  escape  of  atoms  with  a 
thermal  energy  above  the  lattice  potential  well  depth.  We  find  that  the  number  of  atoms 
trapped  in  the  lattice  decays  exponentially,  with  a  time  constant  »  2 (mV GHz) x  A. 
This  scaling  is  consistent  with  loss  caused  predominantly  by  the  escape  of  hot  atoms.  An 
atom  trapped  in  a  potential  well  of  depth  is  heated  at  a  rate  proportional  to  the  photon 

scattering  rate  y, ,  and  will  escape  after  a  time  proportional  to  Upply^E,f  oc  A/F .  Here 

Eg  =  Qikyi2M  is  the  photon  recoil  energy,  and  F  =2;rx5.22MHz  is  the  natural 
linewidth.  Other  loss  processes  that  we  expect  to  contribute  are  optical  pumping  to  the 
F=3  hyperfine  ground  state,  and,  at  long  trapping  times,  escape  in  the  direction 
perpendicular  to  the  lattice  axis. 

An  important  result,  following  fi’om  our  study  of  trapping  dynamics,  is  the 
conclusion  that  resolved-sideband  Raman  cooling  and  quantum  state  control  appears  to 
be  feasible  in  principle.  Sideband  cooling  can  remove  a  quantum  of  vibration  every  few 
oscillation  periods;  this  time  must  be  much  less  than  the  time  required  to  pick  up  a 
quantum  of  vibration  due  to  all  sources  of  heating.  For  the  lattice  parameters  explored  in 
fig.  2,  the  most  rapid  rate  of  increase  in  n  occurs  immediately  after  transfer;  it  is 
approximately  20s‘',  40s''  and  ~  lOOs"'  for  the  3  sets  of  data.  At  the  same  time  the 
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vibrational  oscillation  frequency  in  the  harmonic  approximation  is  42  KHz  and  60  KHz, 
for  lattice 


FIG.  3.  (a)  Typical  time-of-flight  distribution  observed  in  a  simple  time-of-flight 
measurement.  The  corresponding  momentum  distribution  is  indistinguishable 
from  a  Gaussian,  (b)  Time-of-flight  distribution  in  the  presence  of  a  gradient 
magnetic  field.  Each  of  the  9  peaks  correspond  to  atoms  in  a  separate  magnetic 
sublevel. 


depths  of  105  and  209  respectively.  The  increase  in  n  during  an  oscillation 
period  thus  falls  in  the  range  0.5  x  10”^  to  2x10”^,  well  into  the  regime  where  sideband 
cooling  should  be  possible. 
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4,  Stem-Gerlach/time-of-flight  measurement  of  magnetic  sublevel  populations. 

In  order  to  fully  understand  laser  cooling  and  trapping  dynamics  in  near-  and  far- 
off-resonance  optical  lattices,  one  must  characterize  both  the  external  (center-of-mass) 
and  internal  state  of  the  atoms.  To  this  end  we  have  developed  a  new  measurement 
technique,  which  yields  independent  values  for  the  temperature  and  populations  of  atoms 
in  each  ground  state  magnetic  sublevel. 

Our  technique  is  essentially  a  Stem-Gerlach  version  of  the  standard  time-of-flight 
measurement,  which  is  commonly  used  to  obtain  laser  cooled  atomic  momentum 
distributions.  Atoms  are  released  from  the  optical  lattice,  and  fall  to  a  probe  voliune 
located  directly  below  the  lattice;  as  fall  they  experience  a  state  dependent  acceleration 
produced  by  a  magnetic  field  gradient.  The  field  gradient  is  sufficient  strong  to  separate 
time-of-flight  distributions  corresponding  to  atoms  in  different  magnetic  sublevels  (fig. 
3).  A  simple  version  of  this  technique  was  first  explored  in  ref.  [9];  we  have  developed  it 
into  a  precise  measurement  tool  by  (i)  eliminating  precession  of  the  atomic  spin  through 
the  addition  of  a  homogeneous  bias  magnetic  field  at  the  time  atoms  are  released  from  the 
lattice,  and  (ii)  performing  a  detailed  trajectory  analysis  in  the  3-dimensional 
inhomogeneous  magnetic  field. 

Polarization  gradients  of  the  light  field  in  an  optical  lattice  results  in  an 
entanglement  of  atomic  internal  and  external  degrees  of  freedom.  We  have,  therefore, 
proposed  that  the  above  measurement  of  the  atomic  orientation  can  be  used  to  determine 
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atomic  positions  with  sub-wavelength 

measurements  will  be  a  key 
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Fig.  4.  (a)  Populations  (solid  circles)  and  (open  triangles)  in  the  states 
\m  =  ±4/ ,  as  function  of  the  differential  Zeeman  shift  between  those  states  in  units 
of  Upp .  (b)  Temperature  of  atoms  in  the  states  Im  =  4)  (solid  circles)  and  \m  =  -4) 
(open  triangles).  The  data  was  measured  for  a  lattice  with  Upp  =  \A0En  and 
A  =  -20  r . 


ingredient  in  a  series  of  studies  of  macroscopic  quantum  tunneling  in  optical  lattices, 
which  we  are  planning  in  the  near  fuhire. 

5.  Optimized  state  preparation  in  a  near-resonance  optical  lattices. 

Our  research  clearly  demonstrates  that  a  far-off-resonance  optical  lattice  is  best 
loaded  with  atoms  by  transfer  fi’om  a  superimposed  near-resonance  lattice.  In  this 
context,  it  becomes  crucial  to  optimize  the  cooling  process  in  the  near-resonance  lattice, 
so  as  to  maximize  the  occupation  of  the  desired  magnetic-  and  vibrational  ground  state. 
We  have  performed  a  detailed  study  of  the  possibilities  to  enhance  this  initial  state 
preparation,  through  the  addition  of  weak  magnetic  fields.  Our  work  shows  directly  how 
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a  magnetic  field  parallel  to  the  lattice  quantization  axis  modifies  both  the  kinetic 
temperature  and  the  distribution  of  population  among  magnetic  sublevels  (fig.  4). 
Maximum  population  in  one  of  the  two  stretched  states  \m  =  ±4)  is  found  to  occur  when 
the  differential  Zeeman  shift  between  those  states  is  roughly  equal  to  the  modulation 
depth  Upp  of  the  corresponding  lattice  potentials.  We  find  a  maximum  stretched  state 
population  in  the  range  0.70-0.75,  achieved  for  lattice  depths  in  the  range 
Upp  =  70-300  and  a  detuning  A  =  -20  T . 

Because  atoms  are  trapped  near  the  zero-point  of  motion,  the  degree  of  thermal 
excitation  in  the  lattice  is  best  described  by  a  vibrational  temperature.  In  the  context  of 
state  preparation  it  is  even  more  appropriate  to  calculate  and  quote  the  corresponding 
vibrational  ground  state  population  in  the  lattice  potential  wells.  Fig.  5  shows  the  fiiaction 
of  atoms  foimd  in  the  vibrational  grormd  states  of  the  potential  wells  corresponding  to  the 


Fig.  5.  Vibrational  ground  state  populations  (solid  circles)  and  (open 
triangles)  in  the  stretched  states  \m  =  ±4) .  Lattice  parameters  are  Upp  =  10  and 


A  =  -20  r . 
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These  are  the  lattice  parameters  that  produce  the  absolute  maximum  values  of  ~0.3  for  the 
populations.  In  the  absence  of  any  magnetic  field,  we  find  the  vibrational  groimd  state 
populations  for  the  two  stretched  states  to  be  equal,  and  to  attain  a  maximum  value  of 
~0.20,  in  excellent  agreement  with  earlier  theoretical  work  [10]. 

A  magnetic  field  applied  transverse  to  the  lattice  quantization  axis  affects  laser 
cooling  in  a  manner  quite  different  fi'om  an  axial  field.  In  general,  we  find  that  the 
temperature  of  atoms  in  the  stretched  states  shows  a  clear  minimum  for  a  non-zero 
transverse  field;  this  minimum  becomes  even  more  pronounced  in  the  presence  of  an 
axial  field.  The  observed  decrease  in  temperature  is  accompanied  by  a  redistribution  of 
population  among  magnetic  sublevels,  so  that  the  net  result  is  a  decrease  in  the 
vibrational  groxmd  state  populations  of  the  lattice.  We  note  that  the  existence  of  a  global 
minimum  in  temperature  for  B  9^  0  calls  for  caution  if  one  attempts  to  cancel  out 
backgroimd  magnetic  fields  by  minimizing  the  atomic  temperature,  as  is  commonly  done 
in  laser  cooling  experiments. 

6.  Resolved-sideband  Raman  cooling  and  quantum  state  selection. 

We  have  initiated  work  on  resolved-sideband  Raman  cooling,  which  we  plan  to 
implement  using  Raman  transitions  between  magnetic  sublevels,  in  a  manner  closely 
analogous  to  that  demonstrated  for  single  trapped  ions  [2].  To  gain  experience  with 
Raman  spectroscopy  in  a  far-off-resonance  lattice,  we  have  explored  a  simple  ID  version 
of  such  a  scheme  (fig.  6)  [11].  A  pair  of  laser  beams,  of  opposite  circular  polarization 
and  propagating  along  the  lattice  axis,  are  tuned  to  the  "red”  Raman  sideband 


84 


Fig.  6.  Basic  scheme  for  resolved-sideband  Raman  cooling,  (a)  external  Raman 
beams,  or  Raman  coupling  intrinsic  to  the  lattice,  transfers  population  from  state 
to  state  \m=2,n-\).  (b)  Optical  pumping  provides  relaxation  back  to  state 
\m=A,n-\) ,  for  a  net  loss  of  close  to  one  quantum  of  vibrational  excitation. 

optical  pumping.  The  Raman  beams  are  detuned  a  few  linewidths  from  the 
F=4->’F'  =  4  transition,  and  have  adjustable  frequency  difference  5.  This  is 
essentially  the  Raman  cooling  scheme  examined  theoretically  in  ref  [12]. 

We  observe  vibrational  populations  in  the  far-off-resonance  lattice  by  adiabatic 
release  of  the  atoms,  followed  by  a  measmement  of  the  momentmn  distribution. 
Adiabatic  release  maps  the  n  'th  band  onto  the  momentum  interval  («  -l)M  >  p>  nhk 
[13],  and  the  population  of  individual  bands/vibrational  states  can  therefore  be  directly 
measured.  Fig.  7  compares  momentum  distributions,  observed  for  S  detuned  from  or 
resonant  with  the  |m  =  4,  n  =  0/  ->  |m  =  2,  n  =  O)  transition.  It  is  clear  that  atoms  have 
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been  selectively  removed  from  the  \m=A,n  =  Q))  state,  but  have  not  reappeared 
elsewhere.  We  have  not  studied  in  detail  how  atoms  are  lost  during  Raman 
transition/optical  pumping  cycles,  but  it  is  highly  probably  that  optical  pumping  to  the 
F  =3  hyperfine  ground  state  is  responsible. 

In  order  to  minimize  hyperfine  optical  pumping  during  sideband  Raman  cooling, 
it  is  essential  to  minimize  population  in  the  intermediate  excited  state.  This  is  usually 
accomplished  using  far-off-resonance  Raman  beams.  We  have  nearly  completed  a 
theoretical  study  of  a  new  scheme  for  Raman  sideband  cooling,  which  relies  on 
stimulated  Raman  coupling  intrinsic  to  the  lattice  light  field  itself.  A  somewhat 
unexpected  conclusion  is  that  a  destructive  quantum  interference  between  excited 
hyperfine  states  renders  sideband  Raman  cooling  at  best  marginally  feasible,  for  alkali 
atoms  trapped  in  any  1 -dimensional  lattice  configuration.  So  far  this  conclusion  is 
supported  by  our  experimental  efforts,  which  in  the  limit  of  far-off-resonance  Raman 
coupling  have  failed  to  achieve  either  cooling  or  state  selection.  In  a  2D  or  3D  geometry 
the  situation  appears  much  more  promising.  The  extra  degrees  of  freedom  permit  the 
introduction  of  a  7r-polarized  component  of  the  lattice  field,  and  the  resultant  Raman 
coupling  can  be  several  orders  of  magnitude  larger  than  in  a  comparable  ID  far-off- 
resonance  lattice  [11].  We  are  currently  studying  this  situation  in  detail.  We 


Fig.  7.  (a)  Momentum  distributions  measured  after  adiabatic  release  of  atoms 

from  the  far-off-resonance  lattice.  The  solid  line  and  dotted  lines  show  the 
momentum  distribution,  measured  when  the  frequency  difference  S  of  the  Raman 
beams  was  detuned  resp.  resonant  with  the  \m  =4,0)  — » I/w  =  2,0)  transition.  The 
dashed  line  is  the  difference  between  these  distributions,  and  shows  that 
population  has  been  removed  from  the  momentum  interval  -Hk< p<tik, 
corresponding  to  the  ground  vibrational  state,  (b)  Measure  of  the  width  of  the 
momentum  distribution,  as  a  function  of  the  frequency  difference  ^between  the 
Raman  beams.  Arrows  indicate  the  expected  positions  of  the  heating  transitions 
=  4, 0)  Im  =  2,/i  =  {0, 1,2}) ,  and  the  cooling  transition  \m  =  4,l)  ->  |m  =  2,0) . 


have  just  initiated  a  2D  experiment;  demonstrating  resolved-sideband  Raman  cooling  in 
this  configuration  is  a  top  priority,  and  supported  in  part  by  a  new  JSOP  grant  for  the 
period  1997-2000. 

6.  Tunneling  and  quantum  control  in  bipotential  wells. 

A  spin-off  from  our  theoretical  study  of  schemes  for  resolved-sideband  Raman 
cooling,  has  been  the  discovery  of  a  new  and  promising  system  in  which  to  study 
quantum  transport  and  quantum  control.  The  system  consists  of  Cs  atoms  trapped  in  a 
ID  far-off-resonance  lattice,  which  consists  of  a  series  of  double  potential  wells.  Near 
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the  two  minima  of  each  double-well,  the  local  light  polarization  is  predominantly  a+  and 
(T-  respectively,  so  that  bound  states  in  these  wells  have  mainly  m  >  0 ,  respectively 
m  <  0  character.  One  can  then  use  a  small  axial  magnetic  field  to  tune  the  relative  energy 
of  bound  states  in  the  two  wells,  and  thereby  study  both  resonant  and  non-resonant 
macroscopic  quantum  tunneling.  A  measmement  of  the  atomic  magnetization  will  yield 
accurate  information  about  the  atomic  position  at  any  given  time,  and  provide  a  direct 
means  to  observe  tuimeling.  This  will  allow  us  to  address  experimentally  the  effect  of 
dissipation  in  tunneling  [14],  in  a  simple  system  whose  dynamics  can  be  easily  controlled 
and  modeled  via  first-principle  theory.  The  atom/double-well  model  system  has  obvious, 
interesting  parallels  to  other  mesoscopic  quantum  systems,  such  as  quantum  point 
contacts,  SQUID's,  Josephson  junctions  and  small  magnetic  domains. 

8.  Two-dimensional  annular  laser  trap. 

Optical  lattices  are  by  no  means  the  only  configurations  in  which  laser  light  can 
be  used  to  create  traps  for  neutral  atoms.  We  have  developed  a  theory  of  the  two- 
dimensional  motion  of  cold  atoms  in  a  near-resonant  Laguerre-Gaussian  laser  beam  [15]. 
For  a  red-detuned  field  the  laser  beam  provides  an  annular  light  shift  potential  and  the 
atomic  motion  divides  into  vibrational  and  rotational  modes  analogous  to  a  two- 
dimensional  molecule.  In  the  ground  vibrational  state,  this  corresponds  to  a  new 
realization  of  the  two-dimensional  rotator.  We  have  illustrated  the  novel  physics  which 
may  be  explored  with  this  system,  by  showing  that  gravity  acts  analogously  to  a  static 
electric  field  applied  to  a  charged  rotator.  We  have  also  explored  the  possibility  of 
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performing  Raman  spectroscopy  of  the  bound  states  of  the  trap.  Although  it  will  pose  a 
substantial  experimental  challenge,  this  spectroscopy  is  of  considerable  interest  as  it  may 
permit  the  observation  of  new  selection  rules  involving  conservation  of  the  total  orbital 
angular  momentum  of  the  light  and  atoms,  and  may  provide  new  insight  into  the  concept 
of  orbital  angular  momentum  of  light  fields.  We  hope  to  extend  this  work  also  to  laser 
traps  formed  by  e.  g.  Bessel  beams,  in  which  case  the  trapping  potential  has  cylindrical 
symmetry  and  is  approximately  periodic  in  the  radial  coordinate.  The  equivalent 
refractive-index  structures  for  light  waves  have  recently  attracted  great  interest  as 
distributed  feedback  resonators  for  microdisk  semiconductor  lasers. 
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SODIUM  GUIDE  STAR  LASER 
Richard  C.  Powell 
Research  Findings 
Objective 

The  objective  of  this  research  is  to  develop  a  novel  type  of  optically  pumped  solid 
state  laser  system  for  use  as  a  sodium  guide  star  laser  in  adaptive  optics  applications. 

Status  of  Effort 

This  is  a  new  project  that  was  initiated  three  months  ago.  Diuing  this  time,  a 
system  design  concept  has  been  developed  and  materials  and  components  ordered  to 
build  a  laboratory  prototype. 

Accomplishments 

Adaptive  optics  is  an  important  technology  for  military  applications  since  it  can 
provide  diffraction-limited  imaging  from  ground-based  telescopes  for  satellite 
reconnaissance  and  diffraction-limited  laser  beam  delivery  telescopes  for  groundbased 
laser  weapons.  A  “guide  star  laser”  is  one  important  component  of  adaptive  optics 
technology.  So  far,  the  ideal  laser  system  for  generating  sodium  guide  stars  has  not  been 
developed.  The  research  proposed  here  has  the  potential  for  developing  an  all  solid  state 
laser  system  with  the  operating  characteristics  required  for  sodium  guide  star 
applications. 

The  unique  laser  design  proposed  here  utilizes  nonlinear  optics  effects,  such  as 
intracavity  solid  state  stimulated  Raman  shifting  combined  with  second  harmonic 
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generation  to  attain  the  appropriate  operating  wavelength.  The  special  cavity  design 
makes  use  of  the  four-wave  mixing  properties  of  stimulated  Raman  shifting  to  provide 
automatic  quality  enhancement  of  the  beam.  Optimized  high-power  pump  and  Raman 
laser  resonator  designs  will  be  developed  through  computer  modeling  and  prototype 
testing.  The  laser  will  be  capable  of  efficiently  producing  the  pulse  and  energy  format 
required  for  a  sodium  guide  star  laser,  or  for  a  master  oscillator  to  pump  a  power 
amplifier  in  order  to  attain  the  power  levels  required  for  guide  star  applications.  Power 
scaling  remains  the  major  issue.  This  research  will  determine  whether  this  laser  system 
can  operate  by  itself  as  a  high-power  guide  star  laser,  or  if  it  is  better  to  use  it  as  a  seed 
laser  to  pump  high-power  amplifiers. 

The  baseline  approach  we  propose  is  to  fi’equency  double  the  output  of  a  high- 
power  intracavity  solid  state  Raman  laser  already  developed  in  our  laboratory  for  eye-safe 
lidar  applications.  Some  innovative  aspects  of  this  technique  include: 


•  Pump  to  1st  Stokes  conversion  efficiency  ~  90%; 

•  Dual  crystal  second  harmonic  generation  (SHG)  technique  yields  efficiencies  of  ~ 
85%; 

•  Overall  optical-to-optical  efficiency  ~  77%; 

•  Diode  pumping  allows  for  excellent  wall-plug  efficiency  and  minimal  cooling 
requirement; 

•  Compact,  all-solid-state  design  is  amenable  to  rugged  environments; 

•  Excellent  1st  Stokes  beam  quality  (1.1  XDL)  is  obtained  fi-om  a  large  mode- volume, 
highly  multimode  pump  cavity  designs,  which  are  optimized  for  efficient  energy 
extraction,  due  to  intracavity  Raman  beam  cleanup; 

•  "Cleaned  up"  1st  Raman  laser  output  yields  excellent  SHG  efficiency  and  second 
harmonic  beam  quality; 

•  SSRS  has  no  angular,  thermal  of  wavelength  acceptance  products  so  alignment  is 
trivial  and  very  large  (~  5  cm)  crystals  can  be  utilized  to  optimize  the  efficiency; 

•  Intracavity  SSRS  suppresses  laser  "spiking"  in  the  long-pulse  operation; 

•  Intracavity  SSRS  serves  as  an  optical  limiter  to  mitigate  optically  induced  damage. 
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The  Raman  laser  system  will  be  designed  according  to  the  layout  shown  in  Fig.  1. 
This  involves  two  coupled  cavities:  the  pump  laser  cavity  and  the  Raman  laser  cavity. 
The  pump  laser  cavity  is  designed  for  maximum  energy  extraction  to  improve  the  overall 
efficiency  of  the  system.  This  is  achieved  by  requiring  the  laser  medium  (slab  or  rod)  to 
be  the  limiting  aperture  of  the  system.  Intrinsic  to  large  Fresnel  number  cavities  is  a  near- 
field,  highly  multimode,  intracavity  spatial  beam  imprint.  However,  due  to  the 
mechanism  of  intracavity  Raman  beam  cleanup,  none  of  this  higher  order  structure  is 
transferred  to  the  desired  Stokes  output. 

The  laser  resonator  design  is  based  on  a  simple  “L”  -  resonator  with  a  waist 
located  in  the  Raman  medium  as  shown  in  Fig.  1.  An  acusto-optic  modulator  (AOM)  is 
included  to  mode-lock  the  intracavity  field.  A  dual  prism  scheme  is  utilized  to  adjust  the 
cavity  length  to  match  the  natural  driving  frequency  of  the  particular  AOM  employed. 

An  intracavity  etalon  is  placed  in  the  pxunp  laser  cavity  to  confine  the  gain 
bandwidth  to  approximately  1.5  Ghz,  By  this  restriction,  we  are  guaranteed  transform 
limited  pulses  of  temporal  width  of  667  ps.  Our  proven  dynamical  computer  models  have 
established  that  333  ps  transform  limited  Raman  pulses  will  ensue  from  these  pump 
pulses  when  an  appropriate  output  coupling  is  selected.  This  assures  an  output  linewidth 
of  3  GHz,  which  matches  the  inhomogeneously  broadened  absorption  bandwidth  of  a  200 
K  sodium  layer.  The  output  bandwidth  can  be  tuned  to  match  the  varying  diurnal  thermal 
cycles  of  the  sodium  layer  by  a  rotation  of  the  intracavity  etalon. 
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The  central  output  frequency  of  the  pump  laser,  and  hence  system  output,  can  be 
tuned  to  the  D2  resonance  by  injection  seeding  with  a  rapid  tuned,  commercially 
available,  monolithic  isolated  single-mode,  end-pumped  ring  laser  (MISER).  We  propose 
locking  the  output  frequency  to  the  sodium  layer  D2  resonance  by  implementing  a 
positive  feedback  servo-control-loop  to  the  received  backscatter. 


Figure  1.  Raman  sodium  guide  star  laser  system  layout.  The  pump  laser  cavity  and  Raman  laser  cavity  are 
coupled  and  share  a  common  leg  containing  the  Raman  crystal.  The  pump  laser  cavity  is  an  “L”  -  shaped 
cavity  defined  by  two  highly  reflective  (HR)  mirrors  at  the  piunp  wavelengdi.  The  Raman  laser  cavity  is 
mode  matched  by  judiciously  selecting  the  position  and  radius  of  curvature  of  the  end  mirror.  Frequency 
doubling  of  the  Raman  output  radiation  is  performed  in  two  steps  for  an  overall  piunp  to  589  nm  efficiency 
of- 77%. 
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The  Raman  laser  cavity  shares  a  common  leg  with  the  pump  laser  cavity.  The  two 
cavities  are  coupled  by  a  single  dichroic  mirror  that  acts  as  a  high  reflector  (HR)  for  the 
pump  wavelength  and  a  high  transmitter  (HT)  for  the  1st  Stokes  wavelength.  The  mirror 
defining  the  independent  Raman  leg  is  set  to  mode-match  the  two  cavities.  The  first 
fundamental  mode  of  the  Raman  laser  cavity  is  matched  to  the  multimode  beam  radius  of 
the  pump  laser  cavity.  As  a  result  of  intracavity  Raman  beam  cleanup,  the  output  of  the 
Raman  laser  is  guaranteed  to  be  single  mode. 

Our  calculations  show  that  this  system  can  be  scaled  to  tens  of  Watts  with  the 
required  pulse  characteristics  for  sodium  guide  star  applications.  To  reach  a  final  system 
power  of  between  100  W  and  1  kW  will  require  further  research  and  design 
considerations,  and  this  is  one  of  the  major  thrust  areas  of  this  proposal.  Two  possible 
options  to  be  considered  are  using  this  system  as  a  master  oscillator  to  seed  high  power 
amplifiers,  and  developing  new  materials  with  higher  damage  thresholds.  The  current 
technology  of  power  amplifiers  of  the  type  required  for  this  system  have  achieved  powers 
of  over  100  W.  Using  materials  such  as  lithium  borate  for  the  frequency  doubler  will 
significantly  increase  the  optical  damage  threshold  of  the  system. 

There  are  two  basic  Nd  laser  host  -  Raman  laser  crystal  combinations  we  propose 
to  characterize  in  the  feasibility  study.  Our  baseline  approach  is  based  on  Nd:YLF  lasing 
at  the  1.047  mm  transition,  shifted  by  a  KNO3  SSRS  crystal.  We  have  chosen  as  an 
alternative  scheme  a  system  based  on  Nd:YAG  lasing  at  the  1.064  mm  transition,  shifted 
by  a  CaW04  SSRS  crystal.  Each  of  these  schemes  can  produce  a  fundamental  output  of 
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1.1783  mm,  which  when  doubled  is  resonant  with  the  sodium  layer  D2  resonance.  The 
quantum  efficiency  of  each  scheme  is  ~  90%,  so  very  large  conversion  efficiencies  are 
possible.  Both  SSRS  materials  are  commercially  available,  and  have  superior  optical 
quality.  KNO3  has  very  similar  optical  and  SRS  characteristics  to  Ba(N03)2,  which  has 
been  successfully  power  scaled.  CaW04  has  a  relative  low  threshold  for  optical  damage 
and  special  efforts  will  have  to  be  made  to  achieve  the  required  power  scaling  with  this 
material.  We  feel  this  problem  can  be  overcome  through  a  combination  of  material 
improvements  and  system  design  (e.g.  master  oscillator-power  amplifier  configuration). 

Because  of  the  long  pulse  lengths  required  for  efficient  sodium  excitation,  the 
Raman  laser  radiation  has  a  low  peak,  high  average  power.  Efficient  doubling  of  low- 
peak-power  radiation  requires  long  interaction  length,  large  acceptance  parameters,  high 
damage  threshold,  and  low  absorption  in  the  fi’equency  doubling  crystal.  To  achieve 
efficient  doubling  of  low-peak-power  radiation,  a  tight  focus  of  the  pump  radiation  into 
the  nonlinear  crystal  is  necessary  to  obtain  high  peak  intensities.  However,  to  achieve 
doubling  efficiencies  greater  than  50%,  the  bandwidth  of  the  focused  pumps  angular 
spectrum  must  be  narrower  than  the  phase-matching  angular  acceptance  of  the  doubling 
crystal.  This  restricts  the  selection  of  doubling  crystals  and  crystal  cut  geometries.  A 
maximum  angular  acceptance  parameter  and  zero  beam  walk-off  is  achieved  when  phase 
matching  is  possible  in  the  crystallographic  x-y  plane,  i.e.  q  =  90.  In  this  geometry,  the 
phase-matching  becomes  noncritical  (NCPM)  and  a  tightly  focused  beam  can  remain 
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phase-matched  over  very  large  interaction  lengths.  We  are  currently  testing  doped  KTP 
that  possess  the  capacity  to  double  1.1783  mm  to  589.158  nm  in  the  NCPM  geometry. 

Because  all  real  optical  beams  have  a  finite  angular  spectrum,  a  dephasing  occurs 
between  the  fundamental  and  second  harmonic  radiation  that  results  in  a  “back- 
conversion”  of  the  second  harmonic  to  the  fundamental.  For  an  ideal  Gaussian  beam,  this 
occurs  when  60%  of  the  fundamental  have  been  converted  into  the  second  harmonic.  To 
achieve  conversion  efficiencies  greater  than  60%,  a  relative  phase  shift  of  must  be 
imparted  between  the  two  beams  before  further  conversion  can  be  accomplished.  Several 
dual  crystal  -  quadrature  geometries  that  accomplish  a  degree  of  rephasing  have  been 
devised  and  have  been  demonstrated  to  yield  doubling  efficiencies  >  80%. 

The  doubling  technique  we  propose  to  use  utilizes  two  counter-rotated  phase 
plates  between  the  two  crystals  to  rephase  the  fundamental  and  second  harmonic 
radiation.  The  relative  phase  shift  required  to  maximize  the  SHG  efficiency  can  be  tuned 
by  rotating  the  plates.  The  efficiency  is  independent  of  the  pump  drive  and  crystal 
lengths.  Doubling  conversion  efficiencies  of  ~  85%  were  demonstrated  using  this 
technique. 

Scientiflc  Personnel  Supported 

James  Murray 

Publications 


None 
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Interactions/Transitions: 

Collaborations  with  Air  Force  and  Army  Laboratories: 

During  this  funding  period,  we  have  had  site  visit  demonstrations  for; 

Steve  Gottof  and  Cindy  Swim  of  the  Army  Edgewood  Research, 
Development  and  Engineering  Center 
Jay  Fox  of  the  Army  Night  Vision  Laboratory 
Christopher  Clayton  of  the  Air  Force  Phillips  Laboratory 
Ed  Watson  of  the  Air  Force  Wright  Laboratory. 

Each  of  these  individuals  has  expressed  significant  interest  in  the  laser  technology  being 
developed  through  this  project. 

Inventions/Patent  Disclosures 

All  solid  state  guide  star  laser  based  on  Raman  shifting.  (See  attached  form.) 

Honors/Awards 

None 


101 


ULTRA-SHORT  PULSE  PROPAGATION  AND  THE  QUANTUM  STATISTICS 
OF  OPTICAL  SOLITONS 

Ewan  M.  Wright 

Status  of  Effort 

Current  work  in  the  general  areas  of  nonlinear  optics,  in  particular  soliton 
propagation  in  quantum  dot  waveguides,  two-dimensional  motion  of  cold  atoms,  and 
Bose-Einstein  condensation  in  small  atomic  samples  are  progressing  well  with  significant 
recent  progress.  Highlights  include  a  comprehensive  analytical  and  numerical  study  of 
soliton  propagation  in  quantum  dot  waveguides  for  wavelengths  below  the  half  bandgap, 
the  theory  of  artificial  two  dimensional  molecules  formed  fi'om  cold  atoms  trapped  in 
annular  laser  beams,  and  the  prediction  of  collapses  and  revivals  of  the  macroscopic 
wavefimction  for  small  atomic  BECs. 

Accomplishments 

In  the  following  paragraphs  I  will  summarize  recent  accomplishments  in  the  areas 
mentioned  above. 

We  have  recently  been  studying  nonlinear  pulse  propagation  in  quantum  dot 
waveguides  for  wavelengths  below  the  half  bandgap.  The  significance  of  this  work  stems 
firom  recent  interest  in  the  non-resonant  nonlinearity  which  arises  just  below  the  half 
bandgap  in  semiconductors.  These  media,  such  as  InGaAs,  behave  almost  like  Kerr 
media  and  have  allowed  for  the  fabrication  and  demsonstration  of  nonlinear  optical 
switches  (eg.  directional  coupler)  which  operate  near  to  the  ideal  limit,  with  instantaneous 
response  time.  We  have  developed  a  unified  theory  of  pulse  propagation  in  the  vicinity 
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of  a  two-photon  resonance  which  is  applicable  to  semiconductors.  This  theory  accoimts 
for  the  nonlinear  dispersion  in  these  media  which  causes  deviations  from  Kerr-type 
behavior.  For  ultrashort  pulses  the  nonlinear  dispersion  becomes  significant  and  causes 
shock  formation,  which  will  limit  the  nonlinear  optical  switching  capability.  However, 
we  have  found  that  when  combined  with  linear  dispersion  the  nonlinear  dispersion  allows 
for  soliton  propagation.  We  predict  that  these  solitons  should  propagate  in  quantum  dot 
doped  waveguides,  and  that  they  should  allow  for  improved  all-optical  switching  (soliton 
switching). 

In  other  work,  we  have  been  studying  the  2D  motion  of  cold  atoms  in  an  annular 
laser  beam  or  donut  mode.  This  configuration  has  previously  been  studied  but  with  the 
aim  of  trapping  the  atoms  in  the  center  of  the  donut  via  the  dipole  force,  hence  creating  a 
high  density  of  atoms.  In  contrast,  we  have  been  studying  the  case  where  the  atoms  are 
trapped  along  the  peak  of  the  laser  mode,  hence  creating  a  two-dimensional  rotator.  This 
provides  an  atom  optics  realization  of  a  2D  rotator  involving  a  laser  cooled  two-level 
atom  tightly  boxmd  by  the  light  shift  (optical)  potential  supplied  by  a  near-resonant, 
aimular  laser  beam.  This  configuration  is  closely  related  to  previously  demonstrated 
optical  traps  for  laser  cooled  atoms  formed  by  the  light  shift  potential  in  a  tightly  focused 
Gaussian  beam.  Here,  however,  the  light  field  supplies  an  annular  as  opposed  to  a 
Gaussian  light  shift  potential.  Studies  of  atomic  motion  in  periodic  arrays  of  light  shift 
potential  wells  formed  by  laser  standing  waves,  optical  lattices,  have  recently  revealed  a 
close  similarity  to  the  motion  of  electrons  in  a  solid,  and  now  show  promise  as  model 
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systems  in  which  to  study  condensed  matter  phenomena.  Similarly,  the  wavelength  sized 
annular  trap  proposed  here  approximates  a  two-dimensional  molecular  system,  with  a 
spectrum  of  rotational-vibrational  eigenstates.  We  expect  this  "artificial  2D  molecule"  to 
provide  a  new,  well  controlled  model  system  for  molecular  physics,  with  the  obvious 
advantage  that  the  confining  potential  and  spectrum  of  eigenstates  is  exactly  calculable. 
In  its  ground  vibrational  state,  the  system  yields  a  realization  of  the  2D  rotator.  To 
illustrate  the  novel  physics  which  can  be  explored  using  this  system  we  showed  that 
gravity  acts  analogously  to  the  Stark  effect  for  a  charged  rotator  with  an  applied  electric 
field. 

Finally,  we  have  also  been  studying  the  Bose-Einstein  condensation  in  small 
atomic  samples.  In  particular,  we  have  shown  that  the  macroscopic  wavefunction  for 
BECs  composed  of  a  few  thousand  atoms  caimot  act  as  an  order  parameter  for  the  system 
as  it  undergoes  a  series  of  collapses  and  revivals,  reminiscent  of  the  Jaynes-Cummings 
model  of  quantiun  optics.  This  work  has  far  reaching  consequences  for  current 
experiments.  In  particular,  it  calls  into  question  the  role  of  the  macroscopic  wave 
fimction  and  spontaneous  symmetry  breaking  in  determining  the  properties  of  quantum 
gases  such  as  superfluidity.  Even  for  the  large  condensates  now  reported  with  10®  atoms, 
the  collapse  time  is  much  less  than  a  second,  so  that  the  collapses  should  be  observable  in 
current  experiments. 
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Personnel  Supported 

The  following  people  were  supported  by  or  involved  in  JSOP  related  research: 

Ewan  M.  Wright 

Augusto  Rodrigues  (Ph.D  student) 

Keith  Kasunic  (Ph.D  student) 

Kirk  Cook  (Masters  student) 
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Interactions/Transitions 

a.  Conference  Presentations 

The  work  conducted  under  JSOP  has  led  to  several  presentations  at  conferences  which  are 

listed  below: 

1 .  Q.  Feng,  R.  Flesch,  J.  V.  Moloney,  and  E.  M.  Wright,  “Cascade  Breakdown 
Versus  Self-Focusing  for  Femtosecond  Pulses  in  Transparent  Media,”  OSA 
Annual  Meeting/ILS-XI  Program,  Portland,  OR  (Sept.  1995). 

2.  A.  C.  Newell,  Q.  Feng,  J.  V.  Moloney,  and  E.  M.  Wright,  “Eye  Damage:  The 
Role  of  Self-Focusing  in  Laser  Induced  Breakdown  by  Focused  Picosecond 
Pulses  in  Water  and  Other  Materials,”  Nonlinear  Optics  Workshop,  Tucson  (Oct. 
1995). 

3.  M.  Santagiustina  and  E.  M.  Wright,  “Linear  and  Nonlinear  Dispersion  in  Solid- 
State  Solitary-Wave  Lasers,”  contributed  paper,  OSA-IEEE  XI  Advanced  Solid 
State  Lasers  Topical  Meeting,  San  Francisco,  CA  (1996). 

4.  A.  M.  Dunlop,  W.  J.  Firth,  and  E.  M.  Wright,  “Generalized  Mean-Field  or  Master 
Equation  for  Nonlinear  Cavities  with  Transverse  Effects,”  contributed  paper, 
CLEO  96,  Anaheim,  CA  (1996). 

5.  L;  Tomer,  W.  E.  Torraellas,  G.  I.  Stegeman,  and  E.  M.  Wright,  "Trapping  of 
Light  Beams  and  Formation  of  Spatial  Solitary  Waves  in  Quadratic  Nonlinear 
Media,”  invited  paper,  CLEO  96,  Anaheim,  CA  (1996). 


107 


6.  B.  L.  Lawrence,  W.  E.  Torruellas,  G.  I.  Stegeman,  and  E.  M.  Wright,  “Stable  Self- 
Trapping  and  Ring  Formation  in  PTS,”  contributed  paper,  CLEO  96,  Anaheim, 
CA(1996). 

7.  W.  Forysiak,  J.  V.  Moloney,  and  E.  M.  Wright,  “Dynamic  Nonlinear  Optical  Skin 
Effect,”  contributed  paper,  IQEC  XX,  Sydney,  Australia  (July  1996). 

8.  E.  M.  Wright,  D.  F.  Walls,  J.  C.  Garrison,  “Collapses  and  Revivals  of  Bose- 
Einstein  Condensates  Formed  in  Small  Atomic  Samples,”  contributed  paper, 

IQEC  XX,  Sydney,  Australia  (July  1996). 

b.  Consulting 

I  have  been  involved  in  no  consulting  or  advisory  functions. 

c.  Transitions 

As  part  of  my  JSOP  sponsored  work  I  have  been  involved  in  the  development  of 
computer  codes  for  modelling  nonlinear  pulse  propagation  in  water  and  related  media.  In 
particular,  these  codes  have  application  in  studying  laser  eye-damage  mechanisms,  and 
have  been  transferred  to  Brooks  Air  Force  Base  for  that  purpose.  This  work  is  also  part 
of  funding  for  the  Arizona  Center  for  Mathematical  Sciences  through  the  Armstrong 
Laboratory  at  Brooks. 

Discoveries/Inventions/Patent  Disclosures 

None 

Honors/Awards 


None 
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REAL-TIME  ELECTRONIC  HOLOGRAPHIC  MOIRE 
Katherine  Creath 

This  report  is  for  work  performed  between  1 1/1/94  and  6/1/95. 

Objectives 

The  goal  of  this  work  was  to  develop  a  real-time  electronic  holographic  moire 
system  for  simultaneously  measuring  in-plane  and  out-of-plane  displacement  and  stress 
of  mechanical  components. 

Status  of  Effort 

As  reported  last  year,  we  were  able  to  get  some  initial  measurements.  However,  we 
were  limited  by  noise  in  the  measurements  and  were  unable  to  unambiguously  remove 
phase  discontinuities.  Our  subsequent  work  has  concentrated  on  this.  During  the  JSOP 
review  on  2/14/95  I  reported  preliminary  results  showing  how  we  have  been  able  to 
overcome  this  problem.  However,  after  that  review,  the  next  funding  increment  was  not 
approved.  Therefore,  we  finished  up  the  task  we  were  working  on  and  have  not  continued 
this  project. 

Personnel  Supported 

Students  Norbert  Sigrist  and  Joanna  Schmit 


Visiting  Scientist  Peter  Ettl 
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Publications 

The  work  of  the  last  year  is  summarized  in  a  paper  submitted  to  Applied  Optics.  I 
am  currently  making  minor  revisions  asked  for  by  the  reviewers  and  expect  it  to  be 
published  in  about  6  months.  The  title  of  the  paper  is  “Comparison  of  phase  unwrapping 
algorithms  using  gradient  of  first  failure.”  The  reviewers  both  liked  the  idea  of  having  a 
simple  quantitative  means  for  determining  the  success  or  failure  of  a  particular  algorithm. 
There  still  needs  to  be  more  work  done  comparing  algorithms,  but  this  first  step  is  in  the 
right  direction. 

Inventions/Patent  Disclosures 

None 

Honors/Awards 


None 
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DEVELOPMENT  OF  A  GENERALIZED  ASPHERIC  INTERFEROMETER 

John  E.  Greivenkamp 

Objectives 

The  purpose  of  this  program  is  to  develop  an  interferometer  system  that  does  not 
require  use  of  additional  null  optics  to  test  aspheric  optical  surfaces.  Previous  work 
demonstrated  measurement  of  several  hundred  waves  of  departure  from  a  reference 
sphere.  Without  null  optics,  aberrations  are  introduced  at  the  test  part  and  propagate 
through  the  interferometer  (Figure  1),  requiring  calibration  to  relate  the  measured 
wavefront  to  the  part  shape.  This  project  is  an  attempt  to  provide  a  method  of 
interferometric  system  calibration  in  the  presence  of  aspheric  wavefronts.  The  method 
involves  modeling  the  interferometer  using  a  lens  design  program. 

Status  of  Effort 

The  approach  used  is  termed  reverse  optimization.  The  interferometer  optics  are 
measured  as  accurately  as  possible.  This  includes  the  element  spacing  as  well  as  the 
approximate  test  part  location.  This  nominal  interferometer  prescription  is  entered  into 
the  lens  design  program.  The  actual  OPD  measurements  from  the  interferometer  (using 
the  measured  fringe  pattern)  are  entered  as  targets  or  goals  in  the  merit  function.  The 
program  is  then  told  to  optimize  the  system  by  adjusting  the  prescription  to  match  the 
targets.  The  result  of  this  operation  should  give  a  imique  (and  correct)  solution  to  the 
interferometer  design.  The  test  part  shape,  which  is  the  ultimate  imknown,  will  also  be 


computed. 
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Accomplishments 

Aspheric  surface  testing  would  be  greatly  expedited  by  eliminating  the  need  for  a 
null  condition.  To  test  in  a  non-null  configuration,  an  interferometer  must  be  capable  of 
measuring  steep  wavefronts.  These  wavefronts  contain  aberrations  introduced  by  the 
interferometer  optics,  so  non-null  measurements  must  be  calibrated. 

Sub-Nyquist  interferometry  (SNI)  enables  measurement  of  large  wavefronts.  A 
simple  Twyman-Green  interferometer  was  constructed  to  conduct  SNI  and  to  explore  the 
calibration  problem  (Figure  2).  The  characteristics  of  the  interferometer  were 
investigated,  and  wavefronts  with  several  hundred  waves  of  departure  fi’om  a  reference 
sphere  were  recorded  and  successfully  reconstructed.  A  defocused  sphere  was  used  to 
demonstrate  the  need  for  calibration  and  to  serve  as  a  test  subject  for  the  calibration 
study. 

Reverse  optimization  was  used  to  calibrate  wavefiront  measurements. 
Experimental  wavefiront  data  were  entered  into  the  merit  function  of  a  lens  design 
program,  and  optimization  was  employed  to  retrieve  the  prescription  of  the  interferometer 
and  the  part  shape.  The  variables  used  in  the  reverse  optimization  operation  are  shown  in 
Figure  2.  To  simulate  interference  in  the  design  program,  the  reference  wavefiront  was 
stored  in  a  Zemike  phase  surface  on  the  image  plane.  The  multiple  configuration  mode 
of  the  program  was  utilized  to  optimize  several  test  configurations  simultaneously, 
providing  additional  information  to  overcome  a  global  optimization  problem  and  to 
remove  imcertainty  in  the  part  alignment.  Errors  introduced  by  pupil  distortion  were 
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avoided  by  moving  the  stop  to  the  plane  of  the  sensor  in  the  design  program. 
Measurements  from  a  defocused  sphere  were  successfully  calibrated  using  reverse 
optimization.  A  surface  departure  of  1001  was  calibrated  to  better  than  1/4  PV. 

The  results  of  this  calibration  are  shown  in  Figures  3  and  4.  With  no  system 
calibration,  the  measurement  error  is  0.66  waves  peak  to  valley  (Figure  3).  After  the 
reverse  optimization  procedure,  this  error  drops  to  0.12  waves  peak-to-valley  (Figure  4). 

Because  the  interferometer  could  not  be  characterized  sufficiently  for  calibration 
of  measurements  from  aspheres,  simulations  were  conducted.  A  lens  design  program  was 
used  to  develop  an  accurate  model  of  an  interferometer,  including  simulated  fabrication, 
alignment,  and  characterization  errors.  Reverse  optimization  was  applied  using 
wavefronts  generated  by  the  model.  Testing  was  simulated  for  several  conic  asphere 
surfaces.  Reverse  optimization  yielded  calibration  to  1/4  PV  for  conic  aspheres  with 
surface  departures  as  large  as  3001  (Figure  5).  The  results  suggest  that  it  is  feasible  to 
calibrate  non-null  measurements  of  aspheric  surfaces  using  reverse  optimization. 
Personnel  Supported 

Andrew  Lowman  completed  his  Ph.D.  as  a  result  of  this  effort.  The  title  of  his 
dissertation  is  "Calibration  of  a  Non-Null  Interferometer  for  Aspheric  Testing." 
Publications 

1 .  A.E.  Lowman  and  J.E.  Greivenkamp,  “Interferometer  Errors  Due  to  the  Presence 
of  Fringes,”  submitted  to  OS  A  Engineering  &  Laboratory  Notes  (1995) 
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2.  A.E.  Lowman  and  J.E.  Greivenkamp,  “Modeling  an  Interferometer  for  Non-Null 
Testing  of  Aspheres,”  presented  at  Optical  Manufacturing  and  Testing,  Proc.  SPIE 
(San  Diego,  July  1995). 

In  addition,  three  journal  articles  describing  this  work  are  in  preparation  for  submission  to 
Optical  Engineering  and  Applied  Optics. 

Interactions/Transitions 

We  have  had  significant  interaction  with  the  Center  for  Optics  Manufacturing 
(COM)  on  this  project.  COM  is  an  Army-funded  program  headquartered  at  the 
University  of  Rochester.  COM's  primary  funding  has  come  from  the  Army  Mantech 
Program.  This  technology  is  important  to  COM  as  they  develop  automatic  process  for 
the  fabrication  of  aspheric  optics. 

This  work  also  serves  as  the  basis  for  our  proposed  work  under  the  current  ARPA 
Broad  Agency  Announcement  on  Physical  Optics.  That  program  requires  the 
measurement  of  a  large  variety  of  non-rotationally  symmetric  wavefi*ont  and  surfaces, 
and  the  technology  developed  under  this  program  appears  to  be  critical  to  the 
implementation  of  aerodynamically  conformal  windows  to  aircraft  and  missiles. 
Inventions/Patent  Disclosures 
None 

Honors/Awards 


None 
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CONVERGING  LENS 


Non-Null  Test 


CONVERGING  LENS 


ASPHERE 


Figure  1  -  Aberrations  are  introduced  into  the  test  in  non-null  configurations  due  to  the 
change  in  the  path  through  the  optical  system.  In  a  null  configuration  the  rays  retrace  themselves. 
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Figure  2  -  Twyman  Green  interferometer,  with  potential  variable  for  reverse 
optimization  shown.  The  actual  variables  use  are  in  bold. 
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Figure  5  -  Results  of  simulation  showing  the  reverse  optimization  calibration  of 
aspheric  surfaces  with  departures  ranging  fro  25  to  300  waves. 
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WIDE  FIELD-OF-VIEW  MICRO  OPTICS 
T.  D.  Milster 
Research  Findings 
Objective 

This  report  summarizes  the  work  accomplished  for  JSOP  over  the  last  several 
years.  We  started  by  studying  wide  field-of-view  micro-optics.  This  is  an  interesting 
topic,  but  we  found  that  it  is  very  difficult  to  write  submicron  holographic  patterns  on 
curved  substrates.  The  research  effort  then  divided  into  two  areas.  First,  we  investigated 
writing  micro  lithographic  patterns  with  a  near-field  optical  microscope  (Froehlich  et  al, 
1992).  Second,  we  started  to  investigate  other  ways  to  make  lithographic  patterns  on 
curved  substrates,  such  as  using  direct  writing  with  a  focused  laser  beam.  This  part  of  the 
project  eventually  grew  into  a  collaboration  with  Steward  Observatory  for  making 
holographic  null  correctors  for  large-format  telescope  secondary  mirrors  (Vemold  and 
Milster,  1994).  Several  interesting  questions  were  answered,  such  as  how  to  write  a  large 
hologram  without  photoresist  and  how  to  align  optical  and  mechanical  axes  to  submicron 
levels  (Milster  and  Vemold,  1995).  The  major  thrust  of  the  project,  however,  became  our 
research  in  near-field  scanning  optical  microscopy  (NSOM).  The  remainder  of  our  report 
discusses  progress  on  NSOM. 

Status  of  Effort 

Near-field  microscopy  uses  a  subwavelength  aperture  scanned  in  proximity  to  the 
sample  surface.  Light  shines  through  the  aperture  and  illuminates  the  surface  in  a  local 
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region.  The  transmitted  light  is  collected  some  distance  away  from  the  sample  with 
sensitive  detectors.  As  the  aperture  scans  the  sample,  local  variations  in  optical 
properties  produce  modulation  in  the  transmitted  light,  and  a  signal  is  recovered  from  the 
detectors.  An  image  of  the  optical  properties  of  the  surface  is  produced  by  correlating  the 
detector  output  with  the  scan  position. 

Our  NSOM  work  is  divided  into  two  areas,  which  are  modeling  and  experiment. 
The  modeling  effort  uses  a  hybridization  of  finite-difference-time-domain  (FDTD) 
calculations  and  angular  spectrum  calculations.  The  FDTD  code  is  used  to  find  the 
interaction  of  the  light  and  the  sample.  The  angular  spectrum  code  calculates  the  light 
distribution  at  the  detectors.  The  angular  spectrum  calculations  are  straightforward  on  a 
personal  computer,  but  the  FDTD  calculations  require  a  significant  computing  resource. 
The  experimental  work  involves  building  working  NSOM  devices  and  studying  their 
performance.  We  are  currently  on  our  third  generation  device. 

Significant  accomplishments  of  our  modeling  studies  are  listed  below: 
Accomplishments 

1)  BCann  et  al.,  1995a,  show  that  the  output  of  an  NSOM  tip  is  collimated  for  a 
distance  approximately  equal  to  the  aperture  radius.  The  light  then  diverges  linearly.  The 
phase  of  the  near-field  distribution  is  examined  for  the  first  time. 

2)  Kann  et  al.,  1995b,  show  a  detailed  investigation  of  the  relationship  between 
the  sample’s  index  of  refraction,  probe-to-sample  separation,  and  transmitted  power. 
Two  detection  techniques,  which  are  differential  and  total  detection,  are  presented  that 
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each  offer  certain  advantages.  The  differential  technique  has  not  been  described 
previously  and  is  more  robust  to  variations  in  probe-to-sample  separation. 

3)  Kann  et  ai,  1995c,  show  that  a  linear  systems  analysis  of  NSOM  devices  is 
reasonable  over  a  limited  range  of  sample  parameters.  Characteristics  of  the  transfer 
function  are  described. 

4)  Kann  et  al,  1995d,  describe  various  heating  mechanisms  between  the  probe 
and  sample.  By  coupling  a  finite-difference  thermal  model  with  the  FDTD  code,  we 
show  that  the  four  heating  mechanisms  in  the  sample  are:  1)  direct  optical  transfer  and 
absorption;  2)  reradiated  optical  transfer  and  absorption;  3)  thermal  conduction;  and  4) 
thermal  radiation.  For  most  NSOM  probes,  thermal  conduction  and  thermal  radiation 
exhibit  minor  influences  on  sample  heating.  Optical  power  emitted  by  the  probe  and 
absorbed  in  the  sample  plays  the  greatest  role  in  heating  the  sample,  but  reradiated  optical 
power  can  have  a  significant  effect  if  the  polarization  out  of  the  probe  is  oriented  in  the 
plane  of  observation. 

The  significant  accomplishments  of  the  NSOM  experiments  are  listed  below: 

1)  Froehlich  and  Milster,  1994,  show  that  the  minimum  detectable  displacement 

of  an  NSOM  tapered  fiber  probe  using  a  shear-force  servo  is  2.8  x  10'^  A  ^  4Hz  .  This 
is  within  one  order  of  magnitude  of  other  techniques  used  in  high  resolution  atomic  force 
microscopy  (AFM).  Therefore,  the  servo  used  in  om  NSOM  microscope  enables  the 
acquisition  of  both  images  of  optical  surface  properties  with  ultra-high  resolution  and 
images  of  the  surface  topology  that  are  comparable  with  AFM. 
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2)  Froehlich  and  Milster,  1995,  discuss  the  operation  of  the  shear-force  servo 
used  in  our  NSOM  instrument.  Theory  and  experiment  are  in  reasonable  agreement.  The 
tradeoff  between  sensitivity  and  stability  is  described  in  detail. 

3)  In  yet  unpublished  work,  Froehlich  and  Milster  describe  the  nature  of  the  tip 
vibration  amplitude  reduction  that  is  fundamental  to  the  operation  of  shear-force  servos. 
We  first  locked  the  servo  using  the  second  resonance  of  the  probe  tip.  Then  we  added  a 
swept  dither  signal  to  the  probe  tip  near  the  first  resonance  frequency.  The  dither 
vibration  amplitude  versus  frequency  at  the  first  resonance  was  measured.  Results  are 
shown  in  Fig  1(a)  -  (e).  Probe-to-sample  separations  are:  a)  far  from  the  surface;  b)  4.7 
nm;  c)  3.1  nm;  d)  1.9  mn;  and  e)  0.66  nm.  The  increasing  width  of  the  response  curves 
with  closer  spacing  indicates  that  the  reduction  in  probe  dither  amplitude  is  due  to 
increased  damping  and  not  a  change  in  the  resonance  frequency. 

4)  Also  in  yet  unpublished  work,  Froehlich  and  Milster  have  shown  a  dramatic 
example  of  the  usefulness  of  NSOM  when  studying  phase-change  marks  in  optical 
recording  media.  Figure  2(a)  shows  the  topology  of  a  sample  surface  where  600  nm 
diameter  crystalline  marks  are  written  into  an  amorphous  thin-film  layer  on  top  of  a  glass 
substrate.  The  surface  features  around  the  marks  are  due  to  crystallographic  grain 
boimdaries.  The  maximum  deviation  of  the  surface  in  the  mark  regions  is  8  nm,  so  the 
sample  is  essentially  flat.  Figure  2(b)  shows  the  NSOM  image.  Notice  that  no 
correlation  is  observed  with  the  topographical  information  found  in  Fig.  2(a).  The 
striking  optical  contrast  (~50%)  in  the  mark  areas  clearly  shows  the  subsurface  optical 
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characteristics  of  the  thin  film.  In  addition,  the  bright  spot  in  the  center  of  the  marks  is 
not  an  intuitive  result.  Far-field  microscopes,  which  have  much  lower  resolution,  do  not 
show  this  feature.  However,  our  modeling  has  shown  that  the  bright  spot  is  a  real  effect, 
and  its  properties  depend  on  the  shape  of  the  walls  of  the  marks  below  the  surface  of  the 
sample. 

Personnel  Supported: 

Finally,  this  project  has  been  responsible  for  the  graduation  of  two  Ph.D.s  in 
Optical  Sciences.  Dr.  Kaim  now  works  at  USAF  Rome  Laboratory  and  was  sponsored  in 
part  by  the  Senior  Knight  Program.  Mr.  Froehlich  (soon  to  be  Dr.  Froehlich)  will  be 
seeking  employment  after  returning  from  a  bike  tour  in  New  Zealand. 

Publications 

1.  J.  L.  Kann,  T.  D.  Milster,  F.  F.  Froehlich,  R.  W.  Ziolkowski,  and  J.  Judkins, 
“Numerical  analysis  of  a  two-dimensional  near-field  probe,”  Ultramicroscopy  57, 
251-256  (1995a). 
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3.  J.  L.  Kaim,  T.  D.  Milster,  F.  F.  Froehlich,  R.  W.  Ziolkowski,  and  J.  Judkins, 
“Linear  behavior  of  a  near-field  optical  scanning  system,”  J.  Opt.  Soc.  Am.  A 12, 


1677-1682  (1995c). 


124 


4.  F.  F.  Froehlich  and  T.  D.  Milster,  “Minimum  detectable  displacement  in  near- 
field  scanning  optical  microscopy,”  Appl.  Phys.  Lett.  65, 2254-2256  (1994). 

5.  F.  F.  Froehlich  and  T.  D.  Milster,  “Detection  of  probe  dither  motion  in  near-field 
scanning  optical  microscopy,”  accepted  for  publication  in  Appl.  Opt.  (1995). 

6.  T.  D.  Milster  and  C.  L.  Vemold,  “Technique  for  aligning  optical  and  mechanical 
axes  based  on  a  rotating  linear  grating,”  Opt.  Eng.  34, 2840-2845  (1995). 

7.  J.  L,  Kann,  T.  D.  Milster,  F.  F.  Froehlich,  R.  W.  Ziolkowski,  and  J.  Judkins, 
“Heating  mechanisms  in  a  near-field  optical  system,”  submitted  for  publiction  in 
J.  Opt.  Soc.  Am.  A.  (1995d). 

8.  C.  L.  Vemold  and  T.  D.  Milster,  “Non-photolithographic  fabrication  of  large 
computer-generated  diffractive  optical  elements,”  Soc.  Photographic  and  Inst, 
ng.  Proceedings  Volume  2263,  pp.  125-133  (1994). 

9.  F.  Froehilich,  T.  D.  Milster  and  R.  Uber,  “High  resoltuion  lithography  with  a 
near-field  scanning  sub-wavelength  aperture,”  Soc.  Photographic  and  Inst.  Eng. 
Proceedings  Volume  1751,  pp.312-320  (1992). 

Degrees 

J.  L.  Kann,  “Numerical  Modeling  of  a  near-field  optical  system,”  Ph.  D.  University  of 
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F.  F.  Froehlich,  “Resolution  and  contrast  of  optical  and  shear  force  imagery  in  near-field 
scanning  optical  microscopy,”  Ph.  D.  University  of  Arizona,  (final  defense  scheduled 
December  18, 1995). 

Inventions/Patent  Disclosures 
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Figure  Captions: 

1.  NSOM  probe  fi-equency  response  at  the  first  resonance  for  several  probe-to- 
sample  separations:  a)  far  fi'om  the  surface;  b)  4.7  nm;  c)  3.1  nm;  d)  1.9  nm;  and 
e)  0.66  nm.  The  increasing  width  of  the  response  curves  with  closer  spacing 
indicates  that  the  reduction  in  probe  dither  amplitude  is  due  to  increased  damping 
and  not  a  change  in  the  resonance  fi-equency. 

2.  Crystalline  marks  in  an  amorphous  thin-film  background  produce:  (a)  topographic 
scan  showing  ~  8  nm  high  grain  boimdaries  in  the  regions  of  the  600  nm  diameter 
crystalline  marks,  (b)  NSOM  image  showing  the  optical  properties  of  the  marks 
resolved  to  50  nm.  Notice  that  the  topographic  and  optical  images  are  not 
correlated,  which  indicates  that  NSOM  is  producing  new  information  about  the 
surface  optical  properties. 
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Figure  2 


